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Abstract. We consider weight functions Q : C — » R that are locally in a suitable Sobolev space, and 
impose a logarithmic growth condition from below. We use Q as a confining potential in the model 
of one-component plasma (2-dimensional Coulomb gas), and study the configuration of the electron 
cloud as the number n of electrons tends to infinity, while the confining potential is rescaled: we 
use mQ in place of Q and let m tend to infinity as well. We show that if m, n tend to infinity in a 
proportional fashion, with n/m — > f , where < t < +00 is fixed, then the electrons accumulate on 
a compact set St, which we call the droplet. The set St can be obtained as the coincidence set of an 
obstacle problem, if we remove a small set (the shallow points). Moreover, on the droplet St, the 
density of electrons is asymptotically AQ. The growth of the droplets St as t increases is known as 
Laplacian growth. It is well-known that Laplacian growth is unstable. To analyze this feature, we 
introduce the notion of a local droplet, which involves removing part of the obstacle away from the 
set St- The local droplets are no longer uniquely determined by the time parameter t, but at least 
they may be partially ordered. We show that the growth of the local droplets may be terminated in 
a maximal local droplet, or by the droplets' growing to infinity in some direction ("fingering"). 



1. Overview 

1.1. Outline of the paper. In Sections |2] and [3j we study the one-component plasma (Coulomb 
gas ensemble) in two dimensions, and find the quasi-classical limit as the number n of electrons 
tends to infinity while the confining potential is rescaled: mQ replaces Q, where m tends to 
infinity, so that n/m — > t. This was obtained previously by Johansson 1 1 7 ] in the one-dimensional 
context. It turns out that Johansson's proof carries through with only minor modifications also in 
the two-dimensional case, as was explained earlier in our arXiv preprint |13|. Here, we make an 
effort to obtain the result under minimal smoothness and growth assumptions on the potential 

Q. 

In Section 2) we connect the equilibrium measure with an obstacle problem, and show how 
to apply the Kinderlehrer-Stampacchia-Caffarelli theory to obtain a priori smoothness of the 
solutions to the obstacle problem. We also show that the density of the equilibrium measure 
is given by AQ on the droplet, which permits us to reduce the complexity of the equilibrium 
measure to the study of its support (the droplet). The droplet is shown to equal the coincidence 
set for the associated obstacle problem, if we remove the so-called shallow points. For smooth 
strictly convex Q, the topology of the droplets is shown to be simple. 

In Section |5j we introduce the notion of local droplets, which are obtained when we pass 
from the potential Q to its localization Qx for subsets £ c C (cf., e.g. |8|). The local droplets 
are partially ordered, and in Section |6j we study maximal domination chains of local droplets. 
The maximal domination chains either end in a maximal local droplet, or grow to infinity. The 
local droplets appear to be natural from the point of view of physics (see, e.g., ETlO . They are 
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also natural from the mathematical point of view: the description of all possible local droplets is 
exactly the inverse problem of potential theory. 

One purpose with the material on domination chains of droplets in Section [6] is to provide 
a natural setting to analyze Laplacian growth (i.e., the Hele-Shaw equation), which is known 
to be unstable in the forward time direction. This is explained in Section [7] The domination 
chains of droplets are interesting in part because of their integrability nature, especially in the 
case of potentials Q with AQ = constant > near the local droplet (such Q will be called constant 
strength potentials). This will be the topic of a forthcoming paper, where we will discuss the 
algebraic-geometric nature of maximal local droplets for constant strength potentials. 



1 .2. Comments on the exposition. While a few of the results covered in this paper are essentially 
understood, we believe the reader will appreciate a rather self-contained and easily accessible 
exposition. As for the treatment of Johansson's theorems in Section |3j the extension to the two- 
dimensional setting requires some care about details, and as far as we know, no general proof has 
been available so far beyond the arXiv preprint 1 13 1, where an excessive regularity condition was 
made to simplify the presentation (here, we remove that condition by modifying the smoothing 
argument of Johansson's paper [17J; see Subsection l3.2t . 

The connection between equilibrium measures and obstacle problems is known (see, e.g., [20]). 
However, it is perhaps less well known that the Kinderlehrer-Stampacchia-Caffarelli theory (see 
(HI; cf. also 116], where the same technique was used) allows us to develop an understanding 
of the equilibrium measures in terms of their supports, the droplets. This contrasts with the one- 
dimensional theory, where a lot of the difficulty is to determine the density of the equilibrium 
measure. As for the treatment of the Hele-Shaw equation, our approach based on equilibrium 
measures and obstacle problems allows us to develop the theory with low regularity. The standard 
approach to Hele-Shaw flow theory is to use (partial) balayage and variational inequalities, see, 
e.g., [12]. We prefer the obstacle problem approach because it is more intuitive and geometrically 
appealing. 

2. Quasi-classical limit of Coulomb gas ensembles 

2.1. One-component plasma (OCP). In the 2-dimensional Coulomb gas model (or rather OCP, 
the one-component plasma model), we have n electrons located at points {zy}f =1 in the complex 
plane, influenced by an external field. The potential of interaction is 

lo§ ^F' j * k > j ' ke{1 nh 

while the external field potential is denoted by V(z). The function 

V: C -> lU(+oo) 

is lower semi-continuous and sufficiently large to keep the electrons at finite distances. We shall 
supply the precise condition shortly. The combined potential energy resulting from particle 
interaction and the external potential is the function £y : C — » ]R U {+00} given by 

£ v(z) = \ V log ,-^5 + T V( Zj ), z = (z u ..., z„) e C", 

where the summation indices j,k are assumed confined to the set [1, ...,n}. In any reasonable 
gas dynamics model, the low energy states are supposed to be more likely than the high energy 
states. For a positive constant /3, let Z„ = Z n ^y denote the constant 
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where vol 2n denotes the standard volume measure in C" = R 2 ™. We suppose that < Z„ < +oo, 
which means that the potential V imposes a weak localization restraint on the plasma cloud. The 
corresponding Gibbs model then gives the joint density of states 

J_ p -|Sv(z) 

z„ e 

where /3 has the interpretation as the inverse temperature. In terms of the usual van der Monde 
expression 

A(Z) = [I (Zjfc - Zj), 
j,k:j<k 

we may write 

Z„= f |A(z)pe-^/ y ^dvol 2 „(z). 

We thus introduce a probability point process 

Tl n = Un#,v e prob(C") 
(prob(C") is the convex set of all Borel probability measures on C") by setting 



8 v (z) 



dn„(z) = — dvol 2 „(z) = L^Le-iI., vtj) dvol 2n (z), z e C". 

2.2. Marginal measures. For integers k = l,...,n, we define the marginal probability measure 
e protyC*) by setting 

T$\e) = TI n (exC n - k ), 
for Borel measurable subsets e C C k ; in particular, E^"' = FI„. The associated measures 

w _ n\ {k) 
n (n-k)l n 

are known as intensity (or correlation) measures. For k = n, we have T n = n\ Yl n , which is why 
we simplify the notation and write T n = T^. On the other hand, for k = 1, we have (E is the 
expectation operation) 

T*\e) = E[#{; : z y € e}\ 

where it is tacitly assumed that / is confined to the set {1, ... , n], and # denotes counting measure. 
In more explicit form, we have, for n = 2 and k=l, 



(2.1) dlf'CQ = - — — du(Q, C e C, 

where dju(^) = e"2 F ^dvol2(C). More generally, for 1 < k < n and a Borel subset e c C k , we have 

lf\e) = E[#{(ji, ...,j k )e perm(fc, ri) : (z h , ...,z jk )e e}], 
where perm(fc, ri) stands for the collection of all permutations of length k of the set {1, . . . , n}. 
Remark 2.1. In the above definition of the probability measure Tl n , we realize that 

e - j Zj y ( 2 ;)dvol 2 „(z) = djU(Zi) • • • dju(z n ), Z = (Zi, . . . , z„), 

where 

djuCC) = e-5 F(4) dvol 2 (^), «£ e C. 

Most of the above discussion does not depend on this particular structure of the measure (I, 
and we are free to consider more general measures. For instance, this allows us to include the 
one-dimensional theory in the model. 
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2.3. The random normal matrix model. If f> - 2, then the probability measure 



with normalization constant 



dn„(z) = ^ e -WdvoWz) 



Z n = f |A(z)|V^ z <>dvol 2 „(z) 



describes the distribution of the eigenvalues of n x n Random Normal Matrices (RNM) with joint 
probability measure proportional to 

e -bvm m 

where "tr" is the trace, and dM stands for the natural "Haar-type" measure on the submanifold 
of all complex-valued nxn matrices M with M*M = MM* (the normal matrices). In this case the 
point process is determinantal: 

(2.2) dlf (z) = det [K n (z ir Zj )] k i j=i e~ * v ^ dvol 2 ,(z), 
where K„ is the reproducing kernel in the polynomial Bargmann-Fock space 

Pol,, = spanjl, z, . . . , z"' 1 } c L 2 (C, e~ v ). 

We thus consider Pol,, as a finite-dimensional linear subspace of L 2 (C,e~ v ) (linearity is always 
with respect to the field C), and the Gram-Schmidt procedure supplies, for j = 0, . . . , n — 1, 
polynomials pj of degree j and norm 1 such that pj _L pk for j + k. In terms of these orthogonal 
polynomials, we have 

n-1 

(2.3) K n (z, w) = pj(z) pj(w). 



The algebraic mechanism behind the formula for the correlation measure T„ is well understood. 
See, for instance, Mehta's book |19|. 

2.4. Aggregation of quantum droplets. For reasons that will become clearer later on, we shall 
regard the point process T n = n\Tl n (or, equivalently, FI,,) as a quantum droplet. We are interested 
in the transition T n — > T n+ \, which corresponds to adding one more electron to the droplet. A 
direct comparison of the processes T n and r„ + i is not possible, and we are led to consider marginal 
intensities. The following lemma for |3 = 2 has the interpretation that if we add an electron, the 
expected number of A:-tuples of electrons increases everywhere in C*. 



Lemma 2.2. If $ = 2, then 

\/k r (,c) < 



Proof. In view of ||2.3|> , we have 

[k„ +1 (z,-,z ; )]^ =i = [X„(z„z ; )]^ =i + [p n (zi)p n (zj)f. , =i , 

where all matrices involved are positive (semi)definite (the rightmost matrix has rank 1). As we 
compare with (|2.2|>. we realize that the desired assertion 

det[K„(z,-,zy)f . =i < det[x„ +1 (z„z y )J. =i 

is an immediate consequence of the minimax principle (see, e. g., the books of Dunford, Schwarz 
and Gohberg, Krein (TTl). □ 
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Remark 2.3. This "aggregation" property might well be true for all f> < 2 but it certainly fails for 
jS > 2. We consider the illuminating special case H < il / which in the notation of (|2.1t asserts 
that 

(2.4) f || - r/d^Dd^) < 2n(C) f IC-ll^d), C e C. 
Jc 2 Jc 

The measure d/i(|) = e _2y ^'dvoi2(|) can essentially be replaced by an fairly arbitrary positive 
Borel measure (with finite moments). As we plug in the choice df.i = d<5o + dSi, we see that (JZ4) 
is equivalent to 

ICI^ + IC-lf >\, CeC. 

With C - I this gives this gives f> < 2. In fact, it is possible to show that the inequality H 1 ' < 
holds generally for < ft < 2. We outline the argument. It suffices to consider z = in \2A\ . and 
to show that 

(2.5) J [0 + \rjf - || - rtf )d/i(£)d/ifo) > 

for all positive measures [i with finite moments. For < jS < 1 the if triangle inequality shows 
that the integrand on the left hand side is positive point-wise, and the assertion is immediate. 
We turn to the remaining case 1 < p <2. One first establishes with the methods of Calculus that 

(1 + f + 2xf 12 < 1 + t? /2 +fix, - Vf < x < Vt, < t < +co, 
which in complex form becomes 

|1 + zf < 1 + \zf + p Re x, t e D, 
where D denotes the open unit disk in C. By homogenization, this inequality leads to 

II - if < III" + \nf - P min{||p- 2 , \ n t 2 ) Re(77l), |, rj e C, 

so that 

IIP + Irj^ - || - rjf > p min{||p- 2 , \r,t 2 } Re(f]|), |, t ? e C. 
So, to get 02.51 it suffices to obtain 

f mining 2 , 1//" 2 } Re(f?|) d^Od^) > 0. 
Jc 2 

But this is an immediate consequence of Schur 's product theorem for positive definite matrices 
(in this case we have "continuous" matrices), as both min{|||' 3_2 , |r]|^ -2 } and Re(77l) express positive 
definite kernels. 

2.5. Scaling and the class of weights. If we keep the confining potential V fixed, and let n (the 
number of electrons) grow, the process Tl n will generically grow beyond any confinement. For 
this reason, it is necessary to jack up the confinement as n grows. This is achieved by putting 
V - mQ, where m is a scaling parameter and 

Q : C -4RU{+oo) 

is a fixed potential, assumed to be lower semi-continuous. To avoid degeneracy, we must suppose 
that Q < +oo at least on a set of positive area. From well-known physical considerations, it is 
natural to let m be essentially proportional to n. As we are free to pick Q as we like, we may 
assume that the proportionality constant is 1, that is, that m = n + o(n) as n — » +oo. The growth 
requirement on Q which conforms with this normalization is 

(2.6) Q(z) - log |z| 2 -> +oo as |z| -» +oo. 
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2.6. The equilibrium measure. We consider the limit of the point processes 

n 

H m on as n — » +00 while > 1, 

m 

while assuming that Q grows in accordance with | |2.6t . In this case we have convergence of the 
saddle point configurations. More precisely, the probability measures 

(2.7) °* = \L^' 

i 

which minimize the functionals (we write z = (z\, . . . ,z„)) 

# 2 1 v-i 1 2w v-i 

f„, „[ff»] := —, T7 £ «'Q( z ) = T7 > , log : [5 + -; tt > Q(Zj), 

converge as n — > +00 while m - n + o(n) in the weak-star sense of measures to the unique 
probability measure a = Oq which minimizes the weighted logarithmic energy 

(2.8) I Q [o]:= f log-^d t T(5)da(; ] ) + 2 f Qda. 

This comes as no big surprise given the striking similarity of the expressions I* q „[<7«] and Jq[cj]- 
The configuration of points corresponding to a minimizer a n is known as a collection of weighted 
Fekete points, and the measure oq is called the equilibrium measure. The existence and uniqueness 
of the minimizing measure oq is due to Frostman. Let prob c (C) denote the convex body of all 
compactly supported Borel probability measures on C. 

Theorem 2.4 (Frostman). There exists a unique equilibrium measure a = oq such that 

I Q [d] = inf I Q [o], 

a 

the infimum being taken over all compactly supported probability measures 0. 
For the proof, we refer to ||20| . We will write 

(2.9) y(Q):=i Q [ool y*(Q):=y(Q)- f Qdd Q , 

Jc 

for the (modified) Robin constants involved. Let 

(2.10) Lq(& r,) := log + Q(Q + Q(rj), 
and observe that for probability measures a, we have 

(2.11) I Q [a]= f L Q (Ln)d0(i)&0{ii). 

Jc 2 

Next, we introduce the weighted potential 
and observe that since 

(2.12) I Q [a] = £ U^)do(Q 

we expect that the energy minimizer a = oq should have Uq constant on the support 

S = S Q := supp&Q, 

and that constant should also equal the minimum value of IK. We will at times use the notation 
°Q - a lQ] an d Sq = S[Q]. We use q.e. as short-hand for quasi-every where. 



COULOMB GAS ENSEMBLES AND LAPLACIAN GROWTH 



7 



Theorem 2.5 (Frostman). The support Sq of the equilibrium measure oq is compact. Moreover, ify(Q) 
is as in 02.91 1, then U°q > y(Q) cj.e. on C, while U°q < y(Q) at each point o/Sq. The value y(Q) equals 
the minimal energy Iq[&q]. 

For the proof, we refer to ||20l . The number e - ^® is said to be the weighted capacity. 
In terms of the usual logarithmic potential 

we see that for a compactly supported probability measure a, 

U a Q (E,) = f L Q (L n)Ao( n ) = + Q(Q + f Qdo, 
which allows us to write Frostman's Theorem [23] in the following form. 

Theorem 2.6 (Frostman). The support Sq of the equilibrium measure oq is compact. Moreover, ify*(Q) 
is as in $13) . then U d Q + Q> y*(Q) q.e. on C, while U dQ + Q < y*(Q) at each point of Sq. 

Let d mQ,n denote the probability measure o„ given by \2.7\ corresponding to a weighted Fekete 
point configuration (i.e., a minimizing configuration). The convergence to the global energy 
minimizing measure is as follows. 

Theorem 2.7 (Fekete, Totik). We have the convergence 

o m Q,n —* oq as n — > +00 while m = n + o(n) 

in the weak-star sense of measures. Moreover, we have convergence in energy: 

l mQ,n\-°mQ,n] ~> Iq[&q\ = 7(0), as n ~» +°° while m = U + o(n). 

For the proof, we refer to l20l , p. 145. 

2.7. Johansson's marginal measure theorem for the plane. For a probability measure a e 
prob(C) and an integer k = 1, 2, 3, . . ., we denote by a 8i e prob(C ,c ) the product measure given by 

do**(zi, ...,z k ) = da( Zl ) ■ ■ • da(z k ). 

Definition 2.8. We say that Q has extra growth provided that 

(2.13) Q(z) > (1 + So) log(l + |z| 2 ) - Co, ze C, 

holds for some small but positive value of <5o an d some (positive) real constant Co- Moreover, 
we say that Q is regular provided that it is bounded and continuous in an open neighborhood of 

Sq = SUppffQ. 

Theorem 2.9. Suppose Q is regular with extra growth. Then, for every k = 1,2,3, .. ., we have the 
convergence 

H „ dft is n — > +oo while m = n + o(n), 

mQ,n Q x " 

in the weak-star sense of measures. 

Remarkl.lO. (i) Johansson[17JproveshistheoreminthedegeneratereallinecasewhenQ((5) = +oo 
for £ e C \ R (the Hermitian matrix case). This can be viewed as a limit case of our considerations. 
However, the approach of Johansson's proof can be modified so as to include the complex plane 
case stated here. We indicate the necessary modifications in an appendix below, 
(ii) An alternative formulation of Theorem |Z9] runs as follows. As n — > +oo while m - n + o(n), 
the random variables z\, . . . , Zjt on (C", Yl m Q /n ) are asymptotically i.i.d. with law oq. 
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(iii) We now find an application of Theorem |2.9| to linear statistics. Let Cb(C k ) denote the Banach 
space of bounded continuous functions in C k . Moreover, let the trace tr„/ of the function / e Cb(C) 
be given by 

= £ /(*/), 

where the sum as usual runs over j - l,...,n and Z\, . . . , z n are random variables with joint 
probability (C", Tl m Q „). For each j = 1, . . . , n, we have, in view of Johansson's marginal measure 
theorem, for / e Q,(C), 

E[/(z/)] = f f iWn$Q "» f c md& Q (Z)=:(f,& Q ), 
as n — » +00 while m = n + o(n). By forming the average over j, we get, for / e C&(C), 



E 



-tr„/ 
In 



as n — > +00 while m - n + o(n). There is an analogous statement which holds for functions 
/ 6 Cb(C k ) and involves the measure 0^ in place of ctq. This more general statement allows us to 
obtain that for / e Q,(C), (k, k' are fixed integers > 0) 

E 



(^„/) (itr„/) ] ^ {(f,a Q )) k ((f,a Q )f, 



as ;? — » +00 while m = 11 + o(n). Here, as usual, / is the function whose values are complex 
conjugate to those of /. This expresses that ^tr„/ to tends to the constant value (f,OQ) in all 
moments as n — > +00 while m = n + o(n), and hence in particular, we have convergence in 
distribution (as in the weak law of large numbers). 

(iv) We remark that Theorem |Z9] holds independently of the value of the inverse temperature p. 
However, for |3 = 2, much more precise statements have been obtained recently in 0, ||3l . 
The reason why this is possible is the determinantal property l|2.2t . To give some hints about the 
results, we introduce the fluctuation 

f[„f = tr n f -n(f,d Q ). 

In view of (iv), we know that -fl B / — » in moments and hence in distribution as n — > +00 while 
m = n + o(n). Next, suppose n — > +00 while m = n + o(l), which means that m is kept much closer 
to n than before, and suppose also that the function Q is real-analytically smooth with AQ > in 
the interior of Sq (we recall that Sq is the support of the equilibrium measure cjq). In analogy with 
the CLT (central limit theorem), it is shown in UJ, J2] that under some additional assumptions, 
the stochastic variable fl n f converges in distribution to a real-valued Gaussian with expectation 
£f and variance Vf, 

e/ = hj s /a1o§ AQdvo12 ' Vf = hj s |v f |2dvo12 ' 

provided the function f is real-valued, C°°-smooth, and is supported in the interior of Sq. The 
extension to general test functions / is obtained in |3|; the general formulae for ef,Vf include 
boundary effects. 

2.8. Johansson's free energy theorem for the plane. We recall the expression for the normaliza- 
tion constant 



Z m , n = [ lA^M^^dvoWz), 
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which we write in the form 



(2.14) Z m ,„ = f exp | - | £ l q( z ;vZ)c) + |(n - m - 1)£ Q(z y ) idvol 2n (z), 



where Lq is as in 02.101 . The quantity 

1 



■ log Z m , n 



n(n - 1) 

has in the physics literature acquired the name free energy (frequently n 2 is used in place of n(n - 1); 
asymptotically there is no difference). See Definition ^. 8| for the terms regular and extra growth. 

Theorem 2.11. Suppose Q is regular with extra growth. Then 

■ logZ„, /H — > --y{Q) = -^Iq[oq] as n — > +oo ip/zife m = n + o(n). 



n(n-l) ° 4'^' 4 



Remark 2.12. As with Theorem |2.9l Johansson 1 17J proves his theorem in the degenerate real line 
case when Q(£) = +°o for £ e C \ 1R (the Hermitian matrix case). This can be viewed as a limit 
case of our considerations. However, the approach of Johansson's proof can be modified so as 
to include the complex plane case stated here. We indicate the necessary modifications in an 
appendix below. 

2.9. Aggregation of equilibrium measures. We now look at the quasi-classical limit of the 
evolution of quantum droplets (the addition of more electrons to the droplet). This will allow 
us to understand how the quantum process is related to a growth process of Hele-Shaw type for 
compact sets in the plane. 

We restrict our attention to the potentials Q that satisfy a scale invariant version of the growth 
condition j2.6\ , namely 

Q(z) - Alog|z| — > +oo, as |z| — > +oo, 
no matter how big the positive parameter A gets. We will be interested in the evolution of positive 
measures 

ot = QtlQ] := tdQ/t, 

where t ranges over < t < +oo. We write S t = S t [Q] for the support of the measure o>[Q] (i.e., 
St = Sq/j). Note that &t[Q] has total mass t. The process of increasing the parameter t has the 
following interpretation. We consider the limit process of letting n — > +oo while tn = n/t + o(n). 
In other words, m — > +oo while n = mt + o(m). An increase of t therefore has the interpretation of 
increasing the total number of electrons n for fixed m. To rescale, we introduce m' - mt, so that 
the relationship reads n = m' + o(m'). Since mQ = m'Q/t, rescaling also means we must replace 
Q by Q/t. In view of Johansson's theorem, we find that 

#{electrons in e) #{electrons in e) A 
km E = firm E — = to Q/t (e) = a t [Q](e). 

In other words, the growth process dt[Q] is the quasi-classical limit of the growth process of 
adding electrons to the quantum droplet. We shall see that if Q is C 2 -smooth, the measure &t[Q] 
is determined uniquely by its support Sf[Q]. We understand the set St[Q] as a (classical) droplet; 
cf. Subsection 15 . 1 1 for a precise definition. 

Corollary 2.13. The family of measures dt[Q] is monotonically increasing in t. 

Proof. This is true for quantum droplets if jS = 2; the quasi-classical limit does not depend on 
p. □ 

Remark 2.14. It is not hard to write down a potential theoretic proof of this fact; see Proposition 
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3. Appendix: proof of Johansson's marginal measure and free energy theorems for the plane 

3.1 . Fekete configurations. The approach to prove Johansson's theorem in this setting is to show 
that point configurations (zi, . . . ,Z„) whose associated energy functional 

1 ^ 

V ' ),k:j*k 

deviate substantially from the minimum are highly unlikely. We note that since m = n + o(n) is 
assumed, the choice to replace m by n - 1 in the energy is reasonable. Let write 

On '■- &h,(h-1)Q 

for the minimizing (Fekete) measure in the context of Theorem 12. 71 and we also write 

Ini°n\ ■= I n ,(„-X)Q[ n,(n-l)Q] 

for the associated energy. By |20|, pp. 143-145, the sequence of energies l n \d n \ is decreasing in n, 
and converges to /q[6q] = y(Q) as n — » +oo (cf. Theorem|Z7). 

3.2. An entropy estimate. We introduce an auxiliary Borel measurable function <p : C — » [0, +oo) 
with 

(3.1) I (pdvoh = 1, | (Q + | log (p\)(pdvol 2 < +co, 

Jc Jc 

with the understanding that <p log (p = at points where (/> = 0. We sort of artificially smuggle it 
into the expression 11,2.14) for Z mj „: 

Z w ,„ = ( exp L Q (z ; v z^ + ^(n - m - 1) ^ Q(zy) - ^ log 1 [| 0(Z;) dvol 2 „(z). 

Now, by Jensen's inequality, we have, with dcr^ = c/)dvol2, 

(3.2) logZ m ,„ 



- log f { ~ 4 E L q( z /' z ^ + |(» - « - 1)£ Q(Zj) - £ lo g<M z /)} EI ^ dvol 2»( z ) 

C " ^ j,k-j*k j j > j 

f Lq(£, ifjda^da^rf) + ^n(n - m - 1) f Qdcr^ - « f 
Jc 2 2 Jc Jc 



, , ->.— , - ,< r- i - i log^dff^,, 

4 Jc 2 z Jc Jc 

where we used repeatedly that cr^ is a probability measure. We rewrite this as 

1 B n-m-1 C If 

-^logZ^^-^]^^-^- ^Qda,-— ^ logger,, 

This gives (as n —> +oo while m = n + o(n)) 

1 B 
(3.3) hmmf — _ ^ logZ m ,„ > -- J Q [a^J. 

The condition on that cp log <j(> e i. x (C) is of entropy type, and this is the reason why we call i3.3) 
an entropy estimate. We would like to plug in the choice o<p = oq into the entropy estimate < |3.3b to 
get an effective bound. At this point, we do not know enough about oq to be sure whether it is 
of the form a* with <p meeting (13. It . To remedy this, we consider the function <p r : C — » [0, +oo) 
given by (0 < r < 1) 

<M£) = f da Q (r ? ) = f dd Q (4 - t); 
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this amounts to convolution with the normalized characteristic function of the disk D(0, r). The 
corresponding measure 

da,- := do<h r = c/),dvol 2 

is a compactly supported (Borel) probability measure, with density (p r e L°°(C), so that A3. Il l holds 
with cp r in place of <p. By the standard properties of convolutions, a r — > &q in the weak-star sense 
of measures as r — » 0. We claim that we also have convergence in energy, 

(3.4) I Q [a r ] -» I Q [d Q ] = y(Q) as r -> 0. 

Suppose for the moment that we have obtained <|3 .4b - Then we find from the approximation 
procedure that 



(3.5) 



1 B B 

liminf — logZ,„,„ > --IqIoq] = -- y(Q). 



" " 4 « v 4 
To obtain l|3.4t , we note that interchanging the order of integration gives 

Iq[ot]-IqI&q] = 2 f Q(da r -da Q ) + f A r {^ri)da Q {E,)dd Q ^), 
Jc Jc 2 

where 

AM,rf) = -|j f f lo § |(* + T ) - ( ? ? + T ')| - 1{J Sl£ - '7lldvol 2 (T)dvol 2 (T'). 

The support of a r is at most within distance r from the support Sq of 6q, so in view of the 
assumption that Q be bounded and continuous in a fixed neighborhood of Sq, we get 



\ Q(da r - do Q ) = \ Qdo r - \ Qdd Q ^0 
Jc Jc Jc 



(3-6) 

>c Jc 
Next, we rewrite the expression for A, : 



as r — » 0. 



A,(£,rj) 



(T,T')eD(0,r) 2 



log 



1 + 



T - T 



dvol 2 (T)dvol 2 (T') 



= -Ta{ vol 2 (D(0,r)nD(T",r))log 



1 + 



€-1 



dvol 2 (x"). 



We use that the common area of the two intersecting circular disks is 



vol 2 (D(0, r) n D(t", r)) = 2r 2 arccos ^ - r|x"| ^ 1 - 



r"|2 



4r 2 



to get 

ji L r L J T „ 

The identity 



eB(0,2r) 



arccos 



2r 2r 



r»|2 



4r 2 



log 



1 + 



dvol 2 (x"). 



r 

•J-n 



Iog|l + Ae ie |d0 = 27ilog + |A|, A € C, 

where for real x > 0, log" 1 " x = max{0, log x], shows that 

A,.(4,/])=0, if2r<|£-ff|, 

while 

ML n) = ^ L, { arCC ° S £ " Tr "f-^} lo§ j*^jf SdS if 15 - < 2r 
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In the latter case, we may use that for \E, - 1]\ < s <2r, 

2r 



0<arccos--- A /l- — <-, 0<log-— < log 



2r 2r V 4r 2 _ 2' _ & \i - rj\ ~ & |£ - rj\' 
to conclude that 

< A r (£, n) * 8 log n^-j if |5 - r/| < 2r. 
It follows that generally, we have 

2r 



(3.7) 0<A r (£,7 ] )<81og + . 

The measure &q has compact support and finite logarithmic energy, 

so that if we use A3. 7b and the Lebesgue's domintated convergence theorem, we see that 

A r (£, rj) da- Q (5)dff Q (77) -> as r^O. 



f 

Jc 2 



;c 2 

As we combine this with (3.6) . the claimed energy convergence (3.4t is immediate, and hence 
follows. 



3.3. Low probability of high energy configurations. In view of I l3.5t , we have 

(3-8) lo § Zm <" - ~i WQ) + £ )' 

for fixed positive £ and large enough «. In this context, we think of m = m„ as (fixed) sequence 
which depends on n, with m = m n = n + o(n). 
We put 

(1 + |£| 2 )(1 + |7]| 2 ) 

In view of the assumed extra growth 02.131 1, we have 



(3-9) G(£,7/):=logi ', J " > 0, ^eC. 



(3.10) L Q (£, r]) = log + Q(£) + Q(rj) 

> log [ + T^HQtf) + Q(r?)] + log[(l + lif )(1 + |/]| 2 )] - ^ 

14 ~ fjl 1 + °o 1 + oo 



1 + 5 1 + So 

where So and Co are as in d2. 13b - To simplify the notation, we write, with z = (z\, . . . ,z n ), 

»(z) = £ L Q( z /' 2 /c), and G«*»(z) = £ G(z ; ,z,) > 0, 

j,k:j*k j,kj±k 

where it is assumed that j and k range over {1, . . .,n). These expressions are of "double trace 
type" associated with the functions Lq and G (see (12.101 and (13. 9t ). We also have the "trace type" 
expressions (with z = {z\, ... , z„)) 

Q (n) (z) = £ Q(zj) and A<">(z) := £ log(l + |z/). 

i / 

It now follows from (3.101 that 

(3.11) L«(z) > G«*»(z) + ^Q«(z) - -^-n(n - 1), 

« 1 + So 1 + So 
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while a direct application of the extra growth condition | |2,13t leads to 

(3.12) Q<">(z) > (1 + 6 )A <K> (z) - C n. 

The point with introducing this notation is that (|2,14t simplifies to 

(3.13) Z ffl ,„ = £ exp | - ^L«"»(z) + |(n - m - l)Q <H> (z) }dvol 2 „(z), 
while the probability density becomes 

(3.14) dn mQ ,„(z) = -L expj - §L« n »(z) + hn - m - l)Q<">(z)}dvol 2 „(z), 
As mentioned in Subsection 13.11 we have the estimate 



(3.15) -r—r^? & * y<Q)> z = ^ -- z ») eC " 

n(n — 1) y 

We introduce the set 

(3.16) m,e) = [z e C" : -^-L_L«"»(z) < y(Q) + e}, 
where e is a positive real number. 

Proposition 3.1. There exists a positive integer No, which depends one > but not ona>0, such that 

n fflQ ,„(C" \ M(n, e + a))< e-^-W „ > No , 
provided the sequence m— m n = n + o(n) z's kept fixed. 
Proof. By definition, we have 

(3.17) _L^«»»( z )> y( Q) + e + fl/ z6C»\*,e + 4 
n(n — t) « 

We rewrite A3. lit as 

n(n - 1) y (1 + do)n 1 + Oo 

and form a convex combination of 1)3.171 and d3.18t (we keep fixed with < < 1) 

(3.19) — 1— L«"»(z) > (1 - 0)( r (Q) + e + a) + -JL-&q<"> (z) _ 2 C ), zeC"\ Jl(n,e + a). 
n(n — L) w 1 + oq I n J 



The exponent in the density defining Tl m Q,n is (cf. (|3.13t ) 

~\ E L Q( 2 /' 2 *) + f (" - m - !) E Q(z ^ = _ I L Q >>(z) + 5 (n " m " 1 )Q < " > ( 2 )' 
and in view of the estimate (13.191 we get 

(3.20) - ^L«"»(z) + £(n - m - l)Q<">(z) < -^n(n - 1)(1 - 0)(y(Q) + e + a) 

o 

So 



- ^{e(n - 1)^- -(n- m - 1)}q<">(z) + ^n(n - 1), z e C" \ ft(n,e + a). 

If 



(3.21) -^->1- 660 



n-1 1 + §o ' 



14 



HEDENMALM AND MAKAROV 



holds, which is bound to be the case for big enough n (provided 9 is kept away from 0), since 
m = n + o(n), the expression in front of Q^ n \z) on the right hand side of l|3,20t is negative, and we 
may apply | |3.12t to J3.206 , and arrive at 

(3.22) -£ L ^z) + £(n-m-l)Q<»>(z) 

< -tn(n - 1)(1 - 6)(y(Q) + e + a)- ^{d(n - 1)6 - (1 + 6 )(n -m- l))A<">(z) 

+ - 1) - -y-n(n - m - 1), zeC"\ &.{n,e + a). 

As a consequence, we find that 
(3.23) 

n fflQ ,„(C" \ M(n,e + a)) = ^— f exp ( - §L« n »(z) + E( n -m- l)Q< n >(z))dvol 2n (z) 

JC"\m(n,e+a) V 4 ^ J 

< — exp - £n(n - 1)(1 - d)(y(Q) + e + a) + -=£-n(n - 1) - -^n(n - m - 1) 



X | £(1 + |||2)-|(0(»-l)5o-(l+6o)(n- m -l)l dvol2( ^ 

An exercise involving polar coordinates convinces us that for a > 1, 

(3.24) f (1 + m 2 T a dvo\ 2 (Z) = 

Jc « - 1 

and we see that (13.231 entails that 

(3.25) n mQ , n (C"\^(n,e + fl)) 

1 f |S Co/30 C |3 

< — exp - £n(n - 1)(1 - 0)(y(Q) + e + a) + -5£-«(h - 1) - -^-n(n - m - 1) 

In 



B{9(n - 1)6 - (1 + 6o)(« - w - 1)} - 2J ' 



provided that 

m . 06 o 2 



n-1 l + <5 jS(l + 5 )(n - 1) 

Let us assume slightly more, namely that that 

(3.26) ^T>l-^ + 



n-1 l + 6 jS(l + 6 )(n-l)' 

which is a little stronger than J3.2H , and holds for big enough n (as long as 8 is kept away from 
0), since tn = n + o(n). This allows us to get rid of the last factor in the right hand side of J3.25t : 

(3.27) U mQi „{C'\J{{n,e + a)) 

1 { B C Q B6 C Q B ) 

< ^ exp [ - F -n{n - 1)(1 - 0)(y(Q) + e + a) + -^-n(n - 1) - -^n(« - w - 1)}. 

We finally implement the estimate [ 13.81 1, and get 

(3.28) n mQ , n (C n \&(n,e + a)) 

Oy(Q) - (1 - 0)(e + a) + e + 26C - 2C (l "' 



< exp - 1) 



n-1. 
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The constant Q is assumed positive, and we may therefore pick a small 6, < 6 < i, such that 

%(Q) + 2C ]<|. 



Since m = n + o(n), it follows from J3.28t that 
(3.29) "°" 1 ■ 7 " " ' - !/) 



n mQ ,„(C" \ Jl(n,e + a)) < exp \ r -n{n - 1) 



(l-0) fl -(i-0) e + e + o(l))]}. 



Also, by choosing e sufficiently small, we can make sure that 

-(i-e)e + e + o(l)<0 

for big n, so that 03.291 1 gives 

(3.30) n mQ<n (C" \ M(n, e + a)) < exp { - |(1 - 0)an(n - 1) J < exp { - ^an(n - 1) J, 
as claimed. 

3.4. The proof of Johansson's free energy theorem. The claim is that 
1 /S 

(3.31) ^ _ logZ m/ „ -- y{Q) as k > +oo while m = n + o(n). 



Note that by l|3.5t we only need to show that limsup converges to a number < -/3y(Q)/4. To this 
end, we begin by establishing that for < 8 < 1, we have 



(3.32) -- L ^\z) + (n-m-l)Q< n >(z) 



< -- 



■n(n-l)y(Q)- 



(n-1) 



l + <5o 



- (n — m — l) 



Q<»>( z ) + ^„(n-l), 



Q0 



by forming a convex combination of <I3 . 1 5b and (13 . 18b - By applying | |3.12t to d3 . 32b , we get that 
(since the expression in front of Q <n> ( z ) is negative for big m, n with m = n + o{n)) 

(3.33) -lL«"»(z) + (n-m-l)Q<">(z) 

< -l^-nin - l)y(Q) - [(n - 1)5 6 - (1 + 5 )(n -m- 1)](a<">(z) - + ^n(n - 1) 

= -^y%" - l)y(Q) - [(n - 1)6 6 - (1 + 6 )(n - m - l)]A< n> (z) + C dn(n - 1) - C n(n - m - 1). 
We multiply by |S/2 on the left and right hand sides, to get 

(3.34) Z m , n = £ exp { - ^L« n>> (z) + ^(n-m- l)Q<">(z)}dvol 2 „(z) 

< e -§(l-9)«(«-l)y(Q)+|C o 0«(n-l)-Co«(«-m-l)| J" + | < c|2^-|[05 o («-l)-(l+5o)(«-m-l)] c J vo l 2 (^)|" / 

so that in view of (|3.24t , we have 

Z m , n < exp { - £(1 - 6)y(Q) n(n - 1) + ^0n(n - 1) - C n(n - m - 1)}, 
provided l|3.26t is assumed. Taking logarithms, we find that 

^ logZ,„,„ < - J(! - % (Q) + -f 9 - C (l - — ), 
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for big enough m, n with m = n + o(n), since < 13.26t is fulfilled then. As 6, < 6 < 1, can be taken 
as close to as we like, it follows that 

1 j8 
limSUP n(n - 1) lo 8 Zm <" - -jrtQ)' 

The claim is an immediate consequence. □ 

3.5. The proof of Johansson's marginal probability theorem. For a positive real R (a radius), 
we put 

« R (z) = rt{;e{l,...,n}: \ Zj \ > r], 

where (j counts the number of elements, and z = (zj, . . . , z„), as before. We let R be a positive real 
with 

(3.35) 6 log(l + ^)>y(Q) + 2C + l. 

Proposition 3.2. We /zai>e f/ze estimate 

< e, zeM{n,e). 



Proof. We split the integer interval: 

{1, . . . , n] = n(z, Rq) U m(z, R ), n(z, £o) n m(z, R ) = 0, 

where 

n(z,R ) = j/e {l,...,n} : | Z/ | > r], 

so that wr (z) = ji[rt(z,Ro)]- We split the sum defining L*"^(z) accordingly (we use the symmetry 
L Q te,r 1 ) = L Q ( V ,£)): 

L«"»(z)=L I Q (z) + 2Lg(z) + L™(z) 

£ L Q (zj,z k ) + 2 L Q (z p z k )+ L Q (zj,z k ), 

j,kem{z,Ro): j?k j€m(z,R ),k€n(z,R ) jJcen(z,Ro): j±k 

with the obvious interpretation of Lq( z )' ^q( 2 )' anc ^ ^q( z )- From the extra growth condition | [2,13t , 
we see that 

L Q (L rj) > G(4 rj) + 5 log[(l + |£| 2 )(1 + l?]| 2 )] - 2C , 

so that by C335}, 

L Q (£, /]) > <5 log[(l + R 2 )] - 2C > y(Q) + 1, if > Rq. 
This allows us to conclude that 

2Lg(z) + L%(z) > 2n Ro (z)[n - n Ro (z)](y(Q) + 1) + n Ro (z)[n Ro (z) - l](y(Q) + 1), 
As regards the term Lq(z), we may apply (|3.15t to the remaining (n - «R (z))-tuple: 

L l Q {z) >{n- n Ro (z))(n - n Ro (z) - l)y(Q). 
By adding up the terms, we find that 

L«"»(z) = U Q (z) + 2Lg(z) + L™(z) > n(n - l)y(Q) + (n - l)n Ro (z). 
For z e Jl(n, e), we then get 

y (Q) + < _*^f»»( z) < y( Q) + e# 

; K ' n n(n-l) w ; 

from which the assertion is immediate. □ 
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For a point z e C", we define the associated weighted sum of point masses a z e !P C (C) by the 
formula 

(3.36) da z (Q = -Y db z m, £ e C, 

H 

where 6 W means the Dirac point mass at w e C. Also, let Q(C) = C(C) n L°°(C) denote the space 
bounded complex-valued continuous functions on C. 

Proposition 3.3. Suppose o„ = a z is as above, with z = (z\, . . . , z n ) e C". Suppose, moreover, that 

'Q K] = ^1) L<< " >>(Z) y{Q) 
as n — > +oo. TTiefi, as n — > +00, we Tza^e cr„ — > 0q weakly-star. In other words, for each f e Q,(C), we 
have 



I fda„ I /dff Q 
Jc Jc 



as « — > +00. 



Proof. The proof is standard. We choose a weakly-star convergent subsequence, and call the 
limit a*. From the assumptions on the probaility measure a n , we find that almost all its mass is 
concentrated to a fixed compact subset of C (cf. Proposition |3.3l l, and that Iq[o*] < Iq[oq] = y(Q), 
by considering a cut-off of the logarithmic kernel. We leave the details to the interested reader. □ 

Let a> : \l,...,n\ — > be a permutation. For z = {z\,...,z n ) 6 C" we let z m = 

(z w (i),...,z w ( n )) 6 C" be point induced by the permutation. Suppose for the moment that 
/ 6 C;,(C"), and write f a '(z) = f(z a '). By symmetry, we then have 

I / dn,„Q /H = I f° dn m Q /tt , 

JC" JC" 

which gives 

f fdii m Q r n = — I y i / ffl dn m Q„, 

where the sum runs over all permutations co. We next split the integral: 

(3.37) f fdn mQin = - f Vrdn mQ ,„ + i f yrdn mQ/ „. 

By Propositions 13.11 and 13.21 the last term is o(l) as m, n — » +00 while m = n + o(n). In order to the 
understand the remaining term, we should study 

(3.38) on ^ n ' £ )- 

We now focus on the k = 1 case of Johansson's theorem, and restrict our attention to / which 
only depend on the first coordinate, /(z) = /(zi) with some slight abuse of notation. Then (|3.38t 
amounts to the linear statistic 

(3.39) -Y.fbj). {z 1 ,...,z n )ej{{n,e). 

7=1 

By Proposition 13. 2[ only an e proportion of the points Z; may fall outside the disk D(0, Rq), and 
by Proposition [331 the expression J3.39t is close to (the constant!) 

f fd&Q 
Jc 

for small e and large n. The weak-star convergence ^X^q„ oq follows, if we let e approach 
slowly as n — > +00. The remaining case k > 1 is analogous. □ 
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4. An obstacle problem. Smooth potentials 

4.1. Equilibrium measure in terms of an obstacle problem. We consider the cone Sub(C) of all 
subharmonic functions in the plane C, and its convex subset (0 < t < +oo is assumed fixed) 

Sub f (C) := Iv e Sub(C) : limsup [v(z) - Hog|z| 2 ] < +oo|. 

|z|— >+«J ' 

Given Q : C — > 1R U {+°o}, the obstacle problem is to find 

(4.1) Obst,[Q](z) := sup \v{z) : v e Sub f (C) and v<Q on cj. 

Here, we assume of Q - as before - that it is lower semi-continuous, bounded on a set of positive 
area, and that 

(4.2) lim [Q(z)-flog|z| 2 ] = +oo. 

|z|— >+oo 

We think of both Q and t as fixed; we observe, however, that if l|4.2t is fulfilled for one value of f, 
then any smaller positive value works as well. It is easy to check that the supremum in l|4.1t is 
taken over a non-empty collection of functions v (e.g., a large negative constant will satisfy the 
requirements). See e.g. Doob [6| for the potential theory pertaining to obstacle problems of this 
type. For instance, after possibly redefining the function Obst f [Q] on a negligible set (here, this is 
a set of logarithmic capacity 0), we get a subharmonic function. We need to connect the obstacle 
problem ||4jJ with the equilibrium measure theory of Subsections 12. 6l and l2. 91 To this end, let 

ot = o t [Q] := toQ/t, S t - S t [Q] ■= S Q / t = supp at, 

be the scaled equilibrium measure of Subsection |2.9| and its associated support set. We write 

Yt(Q) = ty(Q/t) and f t (Q) = y t (Q) - 1 £ Qdd t . 
For a compactly supported finite positive Borel measure a, let IF denote the logarithmic potential 

and put 

(4.3) Q f (® = y t (Q) - lA-S). 

The function Qt is then subharmonic in C, and harmonic in C \ St, where St = supp 6>. Moreover, 
as it is the total mass of the measure which determines the decay of the logarithmic potential at 
infinity, we have 

(4.4) Q t (z) = t log |z| 2 + O(l) as |z| ^ +oo. 

The following lemma supplies a criterion which allows us to solve the obstacle problem. We 
recall that the logarithmic energy IoM is given by | |2.8t with Q replaced by 0. 

Lemma 4.1. Let a be a compactly supported finite positive Borel measure in C of finite logarithmic energy 
h[o~\ < +°° with total mass \\a\\ = t. Suppose W = c— U° , where c e ]R is a constant. IfW has both 
W < Q q.e. on C and W = Q q.e. on supp a, then W = Obst f [Q] q.e. 

Proof. Without loss of generality, we may assume that Q > 1 on C. The function W is in Sub f (C) 
while W < Q q.e. on C. So W < Q on C \ E, where E c C, is polar (i.e., has logarithmic capacity 
0). Let p be a compactly supported Borel measure on C such that the corresponding potential 
has U p = +oo on E (see, e.g. ||6]|). Put W := (1 - e)W - e'W where e, e' are two small positive 
numbers. If e' is very small (also relative to e), we can make sure that W < Q on C \ E, by using 
that W < Q on C \ E, the standard properties of potentials, and the given properties of Q. Then 
W <Q throughout C automatically as W = -oo on E. We conclude that W < Obst t [Q] on C. By 
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letting e' — » first and then e — » second, we get that W < Obst f [Q] on C \ E, and consequently, 
W < Obstf[Q] q.e. on C. As a side remark, we observe that W is locally of Sobolev class W 1,2 , 
which means that in the sense of distributions, its gradient is locally in L 2 with respect to area 
measure. It remains to show the reverse inequality Obstf[Q] < W q.e. on C. To this end, we 
pick a function v e Sub t (C) with v < Q q.e. on C. It will suffice to show that v < W on C. The 
potential IF is harmonic in C \ S, where S :- supp a, and therefore W is harmonic there as well. 
The assumption on the total mass of a gives that 

(4.5) W(z) = f log|z| 2 + O(l) as |z| +co. 

Next, we consider the difference u = v — W, which is subharmonic in C \ S and has u < q.e. 
on S. Moreover, the assumption that v e Sub t (C) together with 04.51 shows that u is bounded 
from above near infinity. We should like to apply the apply the maximum principle in the open 
set C \ S and obtain that u < on C \ S since u < q.e. on the boundary. However, this is a 
little delicate as the functions are not necessarily continuous up to the boundary. The so-called 
Principle of Domination [20J, p. 104, is a good substitute. To apply it, we need to make the 
technical assumption that v is harmonic in a punctured neighborhood of infinity, because it 
allows us to represent v q.e. in the form b - IT', where b is a constant and v is finite compactly 
supported positive Borel measure. The assumption v e Sub f (C) then gives that v has total mass 
||v|| < t, so that ||v|| < ||a||. From the assumptions we read off that IF < U v + c - b holds q.e. on 
S = supp o, and hence cr-a.e. (because a has finite logarithmic energy), so by the Principle of 
Domination (which again uses that a has finite logarithmic energy, and that ||v|| < ||ct||), we find 
that IF < IF + c - b holds throughout C. The desired conclusion that v < W follows. Next, to 
justify the conclusion that v < W on C holds when we only assume that that v e Subf(C) with 
v < Q q.e. on C, we proceed as follows. If we let e be a small positive real number, and put 

z;(z):=max{^,flog|z| 2 -c}, 

then v e Subf(C) is harmonic in a punctured neighborhood of the point at infinity, and v < Q q.e. 
on C holds if the constant C is big enough positive. So we have the conclusion v < W on C from 
the previous argument. Finally, we first let C — > +oo and afterwards let e — » 0, and obtain v < W, 
as claimed. □ 

We have the following characterization. 
Proposition 4.2. We have 

Obst t [Q](z) = Q t (z), q.e.- z e C. 



Proof. In view of Lemma l4.11 we just need to check that Q t = Q q.e. on St while Qt < Q q.e. on C. 
By Frostman's Theorem |Z6l we have U 5 an + Q/f > y*(Q/t) q.e. on C, while U d Qi> + Q/t = y*(Q/t) 
q.e. on S Q / t . Since d, = td Q/t/ S t = S Q/t , and y*(Q) = ty*(Q/t), this means that IF* + Q > y*(Q) q.e. 
on C, while lF' + Q = y*(Q) q.e. on S t . With Q, = y*(Q) - U d >, this is the same as having Q t < Q 
q.e. on C, while Qt = Q q.e. on St, as needed. □ 

Remark 4.3. The assertion of Proposition ^. 2l is essentially equivalent to that of Theorem 1.4.1 120| . 

We easily recover the density from the potential. 

Corollary 4.4. We have, in the sense of distribution theory, dot = AQ t dA. 
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4.2. The super-coincidence and coincidence sets. We keep the setting of the previous subsection, 
and assume Q:C— >]RU{+oo}is lower semi-continuous, and bounded on a set of positive area, 
subject to the growth condition ||4.2|| . The potential U a is superharmonic for a given finite positive 
compactly supported measure a, and therefore the function Qt defined by l|4.4t is automatically 
subharmonic. In particular, Q t is upper semi-continuous, and we find that the difference Q — Qt 
is lower semi-continuous. It follows that the super-coincidence set 

(4.6) SI = S' t lQ] = jz e C : Q f (z) > Q(z)j 

is compact: it is closed by semi-continuity while l|4.4t and d4.5b show that it is bounded. We note 
that by Proposition 14.21 Qt < Q quasi-everywhere, so that S* t equals, up to a set of logarithmic 
capacity zero, the coincidence set 

jz e C : Q,(z) = Q(z)j. 

In all cases when we have a little regularity, Q t is continuous, and then the super-coincidence set 
S* t is the same as the coincidence set. We therefore refrain from introducing separate notation for 
the coincidence set. 

Proposition 4.5. The function Q t is harmonic in C \ S* t , and as a consequence, St c S* t . In particular, S* t 
is non-empty. 

Proof. We pick a point Zq e C \ S*, so that Q t (z ) < Q(z ). By semi-continuity we get that Q t < Q in 
a neighborhood of Zq. We claim that Qt is harmonic near Zq. If not, we could use Perron's lemma 
and replace Obstf[Q] (which equals Qt q.e., by Proposition 14. 2} on a small disk around zq by the 
harmonic function which has the same boundary values, and get a function which is in Subf(C), 
and bigger than Obst f [Q] while being < Q. This violates the extremality of Obst f [Q], and the 
claim follows. Next, by Corollary l4.41 we see that Zo e C \ St. Since Zq was an arbitrary point in 
C \ S*, the proof is complete. □ 

4.3. A priori smoothness for the obstacle problem for smooth potentials. As before, Q : C — » 
RU j+oo} is lower semi-continuous with l|4.2t where tis a (fixed) positive real. If Q has some degree 
of smoothness, say, e.g., Q : C — » ]R is C 2 -smooth, it is natural to wonder to what extent that carries 
over to Obstf[Q] = Qt. Since we sometimes need to work with slightly less smooth weights, the 
(local) Sobolev classes W 2,p are sometimes more appropriate. By Sobolev imbedding, functions 
in W 2 ' p are continuous provided 1 < p < +oo. The case p = +oo gives the space W 2/0 ° = C 1,1 
of functions whose first order partial derivatives are locally Lipschitz continuous. We use the 
notation W 2 ' p and C 1,1 for local classes unless otherwise stated. 

The following a priori smoothness result is standard in connection with the constrained ob- 
stacle problem discussed below |9], and associated with the names such as Lewy, Stampacchia, 
Brezis, Lions, Kinderlehrer, and Caffarelli. We present the elementary approach recently found 
by Berman JH, which gives the C 1 ' 1 -smoothness part. Berman's approach also applies in the 
several complex variables context. 

Proposition 4.6. IfQe C 2 , then Obst f [Q] 6 C u . More generally, ifQe C u , we still have Obst f [Q] e 
C 1 ' 1 . Finally, ifQe W 2 * for some p,\<p< +co, then Obst t [Q] e W 2 -f . 

Proof. We first show how Q e C 1 ' 1 implies that Obstf[Q] e C 1 ' 1 . We begin by noting that by l|4.4t 
and Proposition H21 

Obst f [Q](z) = Q t (z) < t log(l + |z| 2 ) + C, z e C, 

for a suitable real constant C. By (|4.2|l , the growth of Q(z) is faster than that of Obst t [Q], which 
we can use to show that for some possibly big value of the radius ro, 

(4.7) Obst f [Q](z + w) < Q(z), |z| > r , \w\ < 1. 
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Let 

Mr := sup ||«9 2 {Q(z + tQ}\ : |z| < r, |C| = \, < t < l), 

which is finite for each radius r due to the assumption that Q 6 C ' , and note that by Taylor's 
formula, 

(4.8) Obst f [Q](z + ip) < Q(z + to) < Q(z) + 2 ReMQ(z)] + ±M,.M 2 , |z| < r, \w\ < 1. 
We fix w e C with |zt>| < 1, and put 

Q a ,(z) := iObst f [Q](z + w) + |Obst f [Q](z - a;) - ±M,.>| 2 . 

By a combination of (|4.7[l and l|4.8t , Q w < Q on C, while it is obvious that Q w 6 Sub f (C). So, from 
the definition of the obstacle problem, we see that Q w < Obstf[Q] on C. In other words, 

(4.9) Obstf[Q](z + w) + Obst f [Q](z - w) - 20bst f [Q](z) < M r >| 2 , z e C, \w\ < 1. 
Next, if we divide both sides of < 14,9| | by |w| 2 and then let w — > 0, we get 

5 2 Obst f [Q](z + tC)| f o < M,,, z 6 C, |CI = 1. 

In particular, if z = x + iy, we have 

«9 2 Obst t [Q](z) < M,,, «9 2 Obst f [Q](z) < M, v 

Since Obst f [Q] is subharmonic, that is, 

<9 2 Obst f [Q](z) + <9 2 Obst f [Q](z) > 0. 

holds in the sense of distribution theory, we must then also have 

-M ro < <9 2 Obst t [Q](z) < M,,, -M ro < dpbst t [Q](z) < M,, . 

In particular, then, Obst t [Q] e C 11 . 

As for the remaining case when we have less smoothness, that is, when Q e W 2,p , the asser- 
tion follows from the smoothness theory of constrained obstacle problems (see Lemma 14.71 and 
Theorem |4.9| below). □ 

4.4. A constrained obstacle problem. Let Q : C — » RU j+oo} be lower semi-continuous with l|4.2t 
where t is a (fixed) positive real, as before. Let Q be a (bounded) Jordan domain, and g : dQ. — » IR 
a continuous function with g < Q\aa- Consider the constrained obstacle problem 

ObstQ iP [Q](z) := sup jz?(z) : v e Sub(O), v < Q on O, v = g on <9qJ, zeQ. 

We would like to model the obstacle problem associated with Obst f [Q] in the form of such a 
constrained obstacle problem. The natural way to do this is to put g := Obstf[Q]|^ Q . 

Lemma 4.7. If O is a C^-smooth bounded Jordan domain and g = Obstf[Q]| 5Q , then 

Obst D , e [Q] = ObstrfQ] on Q. 



Proof. We put R = Obst f [Q], R x = Obst,[Q]| Q , and R 2 = Obst Q/P [Q]. The function R a is subhar- 
monic with Ri < Q in Q, and has boundary values Ri\ da — Q- It is now immediate that R\ < R2- 
We proceed to show that R2 < Ri ■ To this end, we let v 6 Sub(O) have D<QonQ and boundary 
data v = g on dQ; we are to check that v < R\. Next, we put v = max{v,Ri}; the function v is in 
Sub(O), has Ri < v < Q on O, and boundary data v\g^ - g. We consider its extension 

J v in Q, 
|R inC\Q. 
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The way things are set up, R < V < Q in C, with V = R on <9Q. We claim that V e Sub(C). It 
is enough to check the mean value inequality along 3D.. For points a e d£l, we have (e > is a 
small real parameter) 

V(a) = Ro(a) < ^- f R (a + ee ie )d6 < f V(a + ee ie )dd. 
2n J-n J-n 

It follows that V e Sub(C) and a fortiori V e Sub t (C) (because of the growth at infinity). We 
conclude that V is a function which we may plug into the optimization problem defining Ro — 
Obstf[Q], and so V < Rq on C. In fact, due to the reverse inequality, we must have V = Ro- In 
particular, v = V\q = R\, and so v < R\. □ 

4.5. Kinderlehrer-Stampacchia-Caffarelli theory. For our purposes it would be enough to con- 
sider the case of C 2 or even C°° potentials Q but since we sometimes have to modify them (see e.g. 
|2 |), the Sobolev classes W 2/P seem to be more appropriate. We generally assume that Q : C — > ]R 
is continuous subject to the growth condition l|4.2} for some (fixed) positive real t. 
We start with a simple observation. 

Lemma 4.8. Let S t = St[Q] and suppose Q e W 2/1 (int St)- Then at is absolutely continuous in int St and 
in fact 

dff f = AQdA on intSf. 



Proof. As we know that do> = AQ t dA in the sense of distributions, and so the same is true if 
we restrict the distributions to the open set int S f , where Q t = Q, and therefore AQ t = AQ as 
distributions. □ 

The following two theorems are adapted from the theory of constrained obstacle problems 
(variational inequalities); this theory is, as mentioned previously, associated with the names of 
Lewy Stampacchia, Brezis, Lions, Kinderlehrer, and Caffarelli, et al. A standard references is l9| , 
Chapter 1 (see also |5|). 

Theorem 4.9. Fix p, 1 < p < +oo, and let D.be a C x '-smooth bounded Jordan domain. We suppose Q 
is W 2 'P -smooth in C, and that g : dQ — > R is a function which is the restriction to dQ of a function in 
W 2 'P(C), with g<Qon dQ. Then Obst Q , e [Q] e W 2 <?(Q). 

Proof. This is explained in Chapter 1 of Friedman's book |9|, see Theorem 1.3.2 and Problem 1 on 
p. 29. □ 

Together with Lemma l4T7l this justifies the W 2 p part of the assertion of Proposition |4.6l Now, in 
view of Proposition ^. 61 if we suppose Q e W 2p for some 1 < p < +oo, the function Q t = Obstf[Q] 
is in W 2 p , and by Corollary 14.41 the measure a> is absolutely continuous (with respect to area), 
and the density is locally in U . By Lemma [4.81 we get 

da, = AQdA on intS t . 

If the boundary dSt has zero area, we can conclude that dot = Is, AQdA. It is remarkable that this 
conclusion holds even when dSt has positive area. 

Theorem 4.10. If, for some 1 < p < +oo, we have Q e W 2 ''', then Q t e W 2/P and 



dot = Is, AQdA. 
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Proof. As S t is the support of at, the measure dt vanishes off St- On St, however, the two W 2,p - 
smooth functions Q and Q t coincide, and by JT81, p. 53, this entails that their partial derivatives 
of order < 2 coincide almost everywhere on St- In particular, AQ< = AQ on St as U functions. In 
view of Corollary l4.41 the assertion is immediate. □ 

Remark 4.11. (a) In the context of Proposition 14.61 it does not help to add more smoothness to Q. 
E.g., if Q e C°° is assumed, we still cannot do better than Obst f [Q] e C 1,1 , at least near dS t . 

(b) The smoothness assumptions of this subsection are excessive, in the sense that it suffices to 
have the required smoothness of Q in a neighborhood of the droplet St. All the statements are 
valid under this weaker assumption. 

(c) By the properties of the 2D Hilbert transform, the assertion that Q t 6 W 2/P is equivalent to the 
property that the density of the absolutely continuous measure do> is in LP (locally). 

4.6. The coincidence set and shallow points. As in the previous subsection, Q : C — » R is 
assumed to be of (local) Sobolev class W 2,p , with 1 < p < +oo. We assume that Q meets the growth 
assumption g2} for all t with < t < T, where T = T(Q) has < T < +oo. 

We recall that we introduced the parameter t to consider the evolution of the renormalized 
equilibrium measures dt = tdQ/t as f moves. The conclusion of Theorem |4.10| allows us reduce the 
complexity and just study the evolution of the droplets S t = S t [Q] = Sq/j. The super-coincidence 
set S* t = S* t [Q] defined by ji.6\ will be referred to as the coincidence set, because the smoothness 
of Q makes Qt continuous. 

We should explain the relationship between the sets St and S* t (we already know that St c S* t ). 
To this end, we say that a point zo 6 S* is Q-shalloiv (with respect to S?) if there exists an open disk 
D centered at Zo such that 

f |AQ|dA = 0. 

Js;nD 

The Q-shallow point in S* form a relatively open subset. We mention in passing that it follows 
from Theorem lilOl that AQ > a.e. on S t . 

Proposition 4.12. The set St is obtained from S* t by removal of all the Q-shallow points. 

Proof. Since Qt and Q are both in C 1,1 and coincide on S* t , we get from |18], p. 53, that AQ f = AQ 
holds a.e. on S*, so that (in the same way as Lemma [4.8| was obtained) 

dd t = l S ;AQdA. 

By comparing with Lemma [4. 81 we see that AQ = a.e. on S* t \ St. To calculate the support of dt, 
we must remove all the points of SJ where there is no |AQ|cL4-mass nearby, that is, the Q-shallow 
points. □ 

4.7. Coincidence sets and the dynamics of droplets. As in the previous subsection, Q : C — > 1R 
is assumed to be of class W 2,p , with 1 < p < +oo. We assume that Q meets the growth assumption 
(|4~2t for all t with < t < T, where T = T(Q) has < T < +co. 

The coincidence set SJ = St[Q] defined by jA.6i is just a little bigger than St (we remove the 
Q-shallow points), but it contains essential information which helps us understand the evolution 
of St as t grows. 

We begin with some elementary properties. 

Lemma 4.13. I/O <h <t 2 <T, then Q h < Q h and, in particular, S* c S* h . 

Proof. Since Sub tl (C) c Sub f2 (C), we clearly have Obst fl [Q] < Obst f2 [Q], from the definition of the 
obstacle problem. The first assertion, Q^ < Qt 2 , now follows from Proposition 14.21 The second 
assertion, the inclusion S* t C S* 2 , is an easy consequence of the first assertion. □ 

Let S* denote the (nonempty compact) set where the global minimum of Q is attained. 
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Proposition 4.14. We have S* c S*Jor allO<t<T. 

Proof. Pick a point a e S* Q , and observe that the function max{Q(a), Q t } is subharmonic (in fact, in 
Sub f (C)) and therefore competes with Q t for the obstacle problem. We conclude that Q(a) < Q f . 
As Qf < Q, it follows that Qt(a) = Q(a), so that a 6 SI. The proof is complete. □ 

Proposition 4.15. IfQ <t\ < t% < T, we have S tl c Sf 2 and < &t 2 - 

Proof. If a point a 6 S * t n SJ 2 is Q-shallow with respect to Sj , then it is also Q-shallow with respect 
to S* ti , since S* t c S* 2 , by Lemma 14.131 The first assertion, Sf a c Sf 2 , now follows from a second 
application of c S* t . The second assertion, < a> 2 , is a consequence of the first assertion 
combined with Theorem 14. 101 □ 

Remark 4.16. This supplies the potential theoretical proof of Corollary l2.13l alluded to in Remark 

EH 

We need the following two lemmas from [20], pp. 227-228. 

Lemma 4.17. (0 < fo < T) The map t i-» S f fs monotonically increasing and left-continuous in the 
Hausdorff metric: 

S t / S to as t / t - 

This means that S to is in a small neighborhood of S f for t < t close to to or, equivalently, that 

S fo = clos^J S f . 

t<to 

Lemma 4.18. (0 < t < T) We ftace 

[I SfCS ?o- 

f <f<r 

In particular, ifS to = S* a , then 

s, = f]s t . 

It is easy to construct examples which show how S* to may contain "seed points" outside the 
main body of St which grow into (small) components of Sf for t > to. The next lemma gives 
a criterion which guarantees that this phenomenon takes place. For a compact set E c C, let 
phull(E) denote its polynomially convex hull, that is, 

phull(E) := |z e C : |p(z)| < max|p| for all polynomials pj. 

The (compact) set phull(£) adds to E all the points of C \ E which belong to bounded connectivity 
components of C \ E (i.e., points invisible to Brownian motion in C \ E starting at oo). 

Lemma 4.19. For all to, t with < to < t < T, w e have the inclusion 

«9[phull(S; o )] c Sf. 

Proof. The standard geometric interpretation of the polynomially convex hull gives that 

5[phuii(s;)]c^s;cs;. 

We need to show that <9[phull(SJ o )] c S t for all t, to < t < T. We argue by contradiction, and 
suppose that there exists a point a e <9[phull(SJ o )] such that a £ C\ St t for some t\, with to < h < T. 
Then a e C \ Sf for all t with to < t < t\, and if t > to is sufficiently close to to, the point a 
belongs to the unbounded component of C \ St. Indeed, choose a small open neighborhood U of a 
avoiding Sfj and since, by assumption, a e <9[phull(SJ Q )], we may assured that there exists a point 
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b e U\ phull(SJ o ). The point b belongs to the unbounded component of C \ S* 1q , so we may connect 
b with oo by a curve y in C \ S* to ; y is at a positive distance from S* tg . By Lemma |4, 181 y C C \ St 
for all t > to close to to, and so the point b - and a fortiori a - is in the unbounded component of 
C \ St- Next, we consider (for t with t > to close to to) the function u = Qt - Qt- Then, by Lemma 
14,131 we have u < 0. Moreover, since a e S* o c SJ, we have Qt {a) = Qt(a) = Q(a), and therefore, 

u(a) = 0. The function Qt is harmonic in C \ St , and, likewise, Qt is harmonic in C \ St, so we 
conclude that u is harmonic in C \ St- The function u then has a local maximum at the interior 
point a, so by the strong maximum principle, we get that u = throughout C \ phull(Sf). This 
does not agree with the known asymptotics ([4.41 . We conclude that the initial assumption must 
be false, so that a e S^ for all t\ with to < t\ < T. □ 

Remark 4.20. The above assertions extend to the case to = if as before S* Q is the set where the 
global minimum of Q is attained, and we put So = 0. 

4.8. Subharmonic potentials. As before, Q : C — » R is assumed to be of class W 2p , so that e.g. 
AQ e i£, c (C). We suppose there exists T = T(Q) with < T < +oo such that 102} holds for 
< t < T while it fails for t > T. 

Lemma 4.21. (0 < t < T) Let D be a bounded domain in C and suppose AQ > in D. Then dD c SJ 
implies D c SJ. 

Proof. The assumption <?D c S* means that Q t = Q on <9D. We write Ro = Qt, and let R\ be the 
function which equals Q in D and equals Qt elsewhere. We observe that Ro < Ri < Q on C, while 
Ro = R-i - Q on 3D. Also, the function Ri is subharmonic. Indeed, ARi = AQ > on D (by 
assumption), and ARi = ARo > on C \ D. It remains to observe that for a 6 dD, 

1 C n 1 f 71 

Ri(fl) = R (fl) < — R (a + ee 10 ) dd < — R^a + ee 10 ) dd, < £ < +co. 

We see that R\ is subharmonic in C, and the conclusion R\ = Ro follows. □ 

Corollary 4.22. IfQ is subharmonic in C, then C\S* t is connected. 

A continuous function h : C — » R is said to be nowhere harmonic if for every open set D c C the 
restriction /z|d fails to be harmonic. 

Corollary 4.23. Suppose Q is subharmonic in C, and that Q is nowhere harmonic. Then S* to c St for all 
t , t with 0<t <t<T. 

Proof. By Corollary l4"!22l the set C \ S* o is connected, and so phull(S* o ) = S* o . By Lemma |4~19l then, 
we arrive at dS* to c St for all t with to < t < T. It remains to check that intSJ o c St for all t with 
t o < t < T. By Proposition 14.121 we just need to show that no point in intSJ o is Q-shallow with 
respect to S* r This is guaranteed by the requirement that Q be nowhere harmonic. □ 

4.9. Convex potentials. We say that a convex function q : C — » 1R is locally uniformly convex if 

l^fA^z) + Re[ed 2 q{z)} > e{z)\E.\ 2 , 5 e C, 

for some continuous e : C — >]0, +oo[. For C 2 -smooth ^, this just says that the Hessian of q is 
(strictly) positive definite everywhere. 

In JT8|, Chapter V, coincidence sets for constrained obstacle problems are considered, and 
under suitable convexity assumptions, the coincidence set is simply connected with C l a -smooth 
boundary (here, < a < 1). The setting is the following. Suppose Q is a strictly convex bounded 
C°°-smooth domain, and let q : Cl — > R be C 2 -smooth and locally uniformly convex, with q > 
on (?Q and mino q < 0. Then, if we put g = in the constrained obstacle problem (see Subsection 
14.411 , tne coincidence set 

S* LXq := jz e Q : Obst Q , [(/](z) = q(z)) 
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is non-empty, compact, simply connected, and equal to the closure of its interior. Moreover, if q 
is C 2,a -smooth for some a, < a < 1, then the boundary dS* Q is a C l,a -smooth Jordan curve, for 
some a', < a' < 1. 

Applied to our setting (cf. Subsection l4.3t , we get Theorem 14. 241 below. Before we formulate 
the theorem, we note that if Q : C \ 1R is convex, and (14. 2t holds for some positive t, then Q must 
grow faster (radially, the growth is at least linear), so that l|4.2t holds for all positive reals t (which 
makes T = T(Q) = +oo). 

Theorem 4.24. (0 < t < T = +oo) Suppose Q : C — » M z's C 2 -smooth and locally uniformly convex with 
l|4.2|l . Then the droplet St is simply connected, and equal to the closure of its interior. Moreover, if Q is 
C 2,a -smooth for some a, < a < 1, then 3S t is a C 1 '"' -smooth Jordan curve, for some a', < a' < 1. 

Proof. We claim that for big enough c, the compact set 

Q c := (z e C : Q t (z) < c] 

is strictly convex with C°°-smooth boundary. In fact, we know from l|4.4t and the fact that Q t has 
the form 

Q t (z) = tlog\z\ 2 + h(z), 

where h is real-valued, bounded, and harmonic in a neighborhood of infinity As c increases the 
sets Cl c cover bigger and bigger portions of the plane C, and the boundary d€l c is contained in a 
fixed neighborhood of infinity for big enough c. The equation defining the boundary is 

\z.\e h(z) = e c , 

and an argument using the harmonic conjugate of /; shows this equation may be written in the 
form 

\z + «o + ci\z~ l + aiz~ 2 + . . . | = e c , 

where the series converges for big |z|. In other words, using the inverse mapping, dCl c is (for big 
c) the image of the circle \z\ = e c under a mapping 

2 I— > z + bo + biz' 1 + fez -2 + . . . , 

which also converges for big |z|. After rescaling by a factor of e~ c , we are talking about the image 
of the unit circle |z| = 1 under the mapping 

zhz + b e~ c + bie~ 2c z~ l + b 2 e~ 3c z~ 2 + 

which for large values of c constitutes a very slight perturbation of the circle |z| = 1, and it is then 
easy to check that the domain inside the curve is strictly convex with C°° -smooth boundary. As 
a consequence, Q c is strictly convex with C°°-smooth boundary for big c. To finish the proof, we 
observe that (cf. Lemma [4.7t 

Obsta c ,o[q] + c = Obst t [Q] on Cl c , 
if q — Q — c. It is immediate that S t = Sq C/9 . The rest follows from Chapter V of IfTHl . □ 

5. Local droplets 

5.1. Localization. We often localize the field Q : C — » ]R U {+00} (which we assume to be lower 
semi-continuous) to a closed set EcC and write 

+00 on <L \ L. 

The function Qz is then also lower semi-continuous. We will assume that Qz meets the growth 
condition l|4.2t , (which is the f-scaled version of {2.6k in case E is compact, this is automatically 
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so irrespective of the behavior of Q near infinity). To avoid triviality, we also need to require that 
< +00 on a set of positive area. We will refer to the closed set E as a localization. 
We will use the notation (which corresponds to the special parameter choice t = 1) 

oq,l = Sf[Q, E] := Qz , Sq iL = S[Q, E] := supp a[Q, E]; 

this conforms with the convention to write &q = d[Q] and Sq = S[Q], We will focus on the 
f-scaled variants (0 < t < +00) 

d,[Q,L] := td Q / tX/ S f [Q,E] := supp<7 f [Q,E] = supp <j[Q/f, E] = S[Q/t,L]. 

We shall also need the modified Robin constant 

y;(Q,Z) = f f (Q z ) 

from Subsection l4.ll 

Lemma 5.1. (0 < t < +00) Suppose E c C is closed, and that Qz meets the growth condition <|4.2|l , while 
Q L < +00 holds on a set of positive area. Then St[Q, E] c E. 

Proof. If the probability measure o~o ■- o[Q/t, E] were to have support outside E, the correspond- 
ing energy 

hz/tloo] = f log 1 da (£,)da (ii) + 7 f Qzda 
Jc 2 1 4 ~ W t Jc 

would necessarily equal +oo, which does not agree with the energy minimizing property of the 

equilibrium measure. □ 

We now compare two different localizations, one contained in the other. 

Lemma 5.2. (0 < t < +00) Suppose Ei, E2 c C are closed with Ei c E2, and that Qj: 2 meets the growth 
condition (|4.2t , while < +00 holds on a set of positive area. IfSt[Q, E2] c Ej, then 

dtlQ, Ex] = o,lQ, E 2 ], S t [Q, Ed = S t [Q, E 2 ], 



Proof. For j = 1,2, we write Oj = d[Q/t, E ; ] and Q ; = Qi^Jt. If Lq is as in l|2.10t , we then get that 
(since Ei c E 2 ) 

L Ql (Z,r}) = L Q2 (Z,ti), (^^eEjxEj, 
and so (since St[Q, E 2 ] = supp a[Q/t, E 2 ] = supp ct 2 c Ei) 

(5.1) Iq> 2 ]= f L Ql (5, ^da^da,^)^ f L Q2 (^,r])dcT 2 (adcT 2 (77)=J Q2 [a 2 ], 

Jc 2 Jc 2 

where we have used the identity 1(2.11)) . As Qi > Q 2 , we have 

where both infima run over a € prob c (C). Combined with [ 15.1| |, this gives Iqj [ffi] < Jqj [a{[, which 
is only possible if o\ = a 2 , by Frostman's theorem (Theorem l2.4t . Finally, if the measures coincide, 
their supports coincide as well. □ 



The typical application of Lemma 15.11 will be when both Ei and E 2 are compact. However, 
already the case when Ei = Sq and E 2 = C is interesting. 

Corollary 5.3. (0 < t < +00) If Q meets the growth condition l|4.2|l , then with E := St[Q] = Sq/j, we 
have 

dt[Q] = d t [Q,n S t [Q,L] = S t [Q]. 
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5.2. Local droplets and the obstacle problem. We let E be a localization, and suppose that 
Q < +00 on a subset of E with positive area. We require that Qz meets the growth condition l|4.2t 
for a positive f, which is kept fixed for the moment (this requirement is void if E C C is compact). 
Then the Borel measure dtlQ, E] is a well-defined positive measure of total mass f, and its support 
Sf[Q, E] is compact with S C E. The following lemma is immediate from Lemma l4~8l 

Lemma 5.4. Suppose Q e W 2 ' 1 (intS) where S = S t [Q, E]. Then cf t [Q, E] is absolutely continuous in 
int S and in fact 

do = AQ dA on int S. 



Proof. This follows from Lemma |4~8l □ 

We are led to the following three definitions. 

Definition 5.5. Suppose Q is in W 2 ' 1 on a neighborhood of S = S f [Q, E]. We say that S is a local 
(Q, t)-droplet with localization E if the following equality holds (a = dt[Q, E]): 

da = l s AQdA. 



Definition 5.6. A compact set S c C is a local (Q, t)-droplet if it is a local (Q, t)-droplet with respect 
to some localization E. 

Definition 5.7. A compact set S c C is a global (Q, t)-droplet if it is a local (Q, f )-droplet with respect 
to the localization E = C. 

Remark 5.8. (a) There is at most one global (Q, t)-droplet S, as it is given by S = St[Q]. Theorem 
14.101 guarantees that it exists under the growth requirement l|4.2t and the additional regularity 
Q e W 2,p (this is true also if Q is in W 2,p only in a neighborhood of S). In contrast, there may exist 
several local (Q, £)-droplets. 

(b) If the boundary dS of the set S = S[Q, E] has zero area, then S is a local (Q, t )-droplet if and 
only if and only if a — dt[Q, E] is absolutely continuous. 

(c) Two local (Q, t)-droplets Si, S2 cannot have the containment Si C S2 unless Si = S2. 

(d) The point with the definition of local droplets is that we may focus on the support S = supp 
rather than the (generally more complicated) equilibrium measure o (with respect to the weight 
Ql). 

(e) In the above definition, it is possible to weaken the smoothness assumption on Q to W A1 
smoothness, which just asks that the function and its Laplacian are both locally integrable. 

Proposition 5.9. If S is a local (Q, t)-droplet with respect to some localization E = Eo, then it is a local 
(Q, t)-droplet with respect to the minimal localization E = S. 

Proof. This follows from Lemma [5 . 2 1 with Ei = S and E2 = Eo. □ 

Remark 5.10. If S = S t [Q, E] is a local (Q, f)-droplet, then the associated measure a = a t [Q, E] is 
absolutely continuous. It is possible that the converse might be true (cf. Lemma [5 .41 above) . For 
the moment, we have a weaker statement. Suppose first that Q is in ]N 2, f in a neighborhood of S, 
for some p, 1 < p < +00. The statement now runs as follows: if a is absolutely continuous with 
density in U for some p, 1 < p < +00, then S is a local (Q, £)-droplet. Indeed, from the properties 
of the 2D Hilbert transform, we get that the function 

(&) t ® = y ;(Q,E)-um 

is in W 2,p and from Proposition ^. 5| we have that 



(QzMQ = Q(Q, Z e S, 
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so that by Il8l , p. 53, we get 

A(&M£) = AQ(Q, I 6 S, 
as distributions, which leads to the desired result. 

For a compact ScC, we define the corresponding (weighted) logarithmic potential 

(5.2) U& s (£) = £log — !-p AQ(n)dA(rj). 

We have the following characterization of local (Q, l)-droplets. We recall the notion of Q- 
shallow points from Subsection |4.6l 

Theorem 5.11. Suppose that S c E c C, where S is compact and E closed, and that Q is in W 2 ' 1 in a 
neighborhood ofS. Then S is a local (Q, t)-droplet with localization E if and only if: 

(i) AQ > a.e. on S, 

(ii) S contains no Q-shallow points, 

(iii) j s AQdA = *' 

(iv) U S, Q + Q = y]{Q, S) a.e. on S,for some real constant y* t (Q, S) (the modified Robin constant), and 

(v) U s -Q + Q>y t (Q,S)q.e.onL. 

Proof. We first establish the necessity of conditions (i)-(iv). So, we suppose that S is a local (Q, f)- 
droplet. As da = lgAQdA is positive with mass t, and S is its support set, conditions (i)-(iii) are 
necessary. The necessity of condition (iv) and (v) follows from Frostman's Theorem 12.61 (with 
Qs/t in place of Q, where S is used as a localization). 

We turn to the sufficiency of the conditions (i)-(v). We write W :- c - U^' S , where the constant 
c = Y]{Q, S) is as in (iv). By (iv), we then have W = Qz q.e. on S while (v) gives W < Qz q.e. on 
C. By LemmaED we get that W = Obst f [Qi;]. Next, Proposition S2] and Corollary S3] show that 

d& t [Q,S] = AWdA = -AU as dA = l s AQdA. 

So we have a local (Q, f)-droplet with localization E. □ 

Remark 5.12. To characterize the local (Q, t)-droplets, we use the minimal localization S. We see 
that condition (v) becomes vacuous and may be removed. 

Corollary 5.13. Suppose that S c E c C, where S is compact and E closed, and that Q is in W 2,1 in a 
neighborhood of S. Consider the function 

Qs:=f t (Q,S)-U& s , 

where y* t (Q, S) is the constant in Theorem \5.11\ Then S is a local (Q, t)-droplet with localization E if and 
only if. 

(i) AQ > a.e. on S, 

(ii) S contains no Q-shallow points, 

(iii) j s AQdA = l ' 

(iv) _Qs = Q q-e- on S, and 

(v) Q s < Q q-e. on E. 

Moreover, if (i)-(v) are assumed, then Q s e Sub t (C) is harmonic onC\S, with asymptotics 

Q s (z) = nog|z| 2 + 0(l) as Izl^+cx,. 

As a consequence, we have q.e. 

Qs = (&j t = ObsttKM 
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Moreover, if, for some p with 1 < p < +00, we have Q 6 W 2,p in a neighborhood of S, then Qs e W 2,p as 
well. 

Proof. It is clear that from the properties of logarithmic potentials that U® ,s is subharmonic in 
C and harmonic in C \ S, with the corresponding asymptotics at infinity as a consequence of 
condition (ii) of Theorem 15.111 Moreover, the properties of the 2D Hilbert transform show that 
if Q e W 2,p in a neighborhood of S, then U®' S e W 2,p , for 1 < p < +00. These properties are then 
inherited by Q$. □ 

If there is some room to wiggle between the set S f [Q, E] and the localization E, then the set 
Sf [Q, E] is automatically a local (Q, £)-droplet: 

Theorem 5.14. Suppose Q e W 2 * 'for some p, 1 < p < +00. If for a localization E we have S = S t [Q, E] c 
int E, then S is a local (Q, t)-droplet with localization E. 

Proof. This is Theorem 14 . 1 1 for Qn in place of Q. □ 
Remark 5.15. 

The modified Robin constant y*(Q, S) may be written out explicitly: 
(5.3) f t (Q,S) = ] f log— L-AQ(£)AQ(r ] )cL4(£)dA(r ] ) + f QAQdA. 

5.3. Characterization of local droplets. We need the concept of local Q-droplets. We consider 
compact localizations E only which means that no requirement on Q near infinity is needed, just 
that Q:C— >]RU{+oo}is lower semi-continuous and has Q < +00 on a subset of E with positive 
area. We recall the concept of a (Q, £)-droplet, which presupposed that Q was W 21 -smooth near 
S. 

Definition 5.16. A compact set S C C is a (local) Q-droplet if it is a local (Q, f)-droplet for some t 
with < t < +00. 

We see that Theorem 15.11 l has the following consequence. 

Corollary 5.17. Suppose that S c C is compact, and that Q is in W 2,1 in a neighborhood ofS. Then S is 
a local Q-droplet if and only if: 

(i) AQ > a.e. on S, 

(ii) S contains no Q-shallow points, and 

(iii) U S, Q + Q is constant q.e. on S. 

By Sobolev imbedding, we have W 2,p C C 1 for 2 < p < +00. The following characterization 
will prove useful later. 

Proposition 5.18. (0 < t < +00) Suppose S c C is compact with S = clos int S, and that Q e W 2,p in a 
neighborhood of S, for some p, 2 < p < +00. We then have: 

(i) IfS is a local Q-droplet, then d(U Q - s + Q) = on S. 

(ii) If S is connected and d(U® ,s + Q) = on dS, and if S has no Q-shallow points, then S is a local 
Q-droplet. 

Proof. We first treat part (i). So, we assume that S is a local Q-droplet. By Corollary 15. 171 U®' S + Q 
is constant q.e. on S. As both Q and !i Q,s are in W 2 ' p in a neighborhood of S, we conclude from 
ED, P- 53, that 5(U Q - S + Q) = a.e. on S. By Sobolev imbedding, d(U Q - s + Q) is continuous in a 
neighborhood of S, and so d( U^ ,s + Q) = on int S and a fortiori (by the topological assumption) 
on S. 

We turn to part (ii). Consider the function F := 5(U® ,S + Q), which is in W 1/P in a neighborhood 
of S, and therefore continuous. We have 

dF = A(!i Q ' S + Q) = -l s AQ + AQ = a.e. on S. 
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Hence F is conjugate holomorphic in the interior of S and since F = on the boundary, we have 
F = on S. If S is connected, then this implies that Us + Q is constant on S, so by Corollary l5.171 
S is a local Q-droplet. □ 

6. Chains of local droplets 

6.1. A partial ordering of local droplets. We recall that S is a local Q-droplet if it is a local 
(Q, £)-droplet for some t with < f < +oo. For the concept to make sense, we need to ask that 
Q : C — » ]R U {+00} is lower semi-continuous and W 2,1 -smooth near S. Given a local Q-droplet S, 
the corresponding value of (the evolution parameter) t is easily calculated: 

We note that by Corollary |5.171 AQ > on S. To simplify the presentation, we shall assume that Q 
is W 2 ' p -smooth in Cfor some p, 1 < p < +oo. 

Lemma 6.1. Let S2 be a local Q-droplet, with ti = t(Q, S 2 )- Ifh has < t\ < t 2 we put Si := S tl [Q, S{\- 
Then Si is a local Q-droplet, with t\ = t(Q, Si). 

Proof. We should study the measure o~\ := a tl [Q, S2], which by Proposition [42] and Corollary 14.41 
is obtained from ObsttJQsJ by applying the Laplacian. From t\ < tz and the definition of the 
obstacle problem, we see that 

Obst tl [Q S2 ]<Obst t2 [Q S2 ] = Qs 2 , 

where we use Corollary 15 . 1 31 to get the rightmost identity. A moments reflection, using that 
Qs 2 < Q, reveals that in fact 

Obst fl [Q S2 ] = Obst t2 [Q s J. 

Since AQ is in If locally, Qs 2 is W 2,p -smooth, and by Theorem 14 . 1 1 with Qs 2 in place of Q, we get 
that Si is a local Q-droplet. □ 

Lemma [6.11 allows us to introduce a partial ordering in the set of all local Q-droplets. 

Definition 6.2. Let Si, S2 be two local Q-droplets, and write fy = f(Q, Sy), j = 1, 2. We write Si < S2 
if Si c S 2 andSi = S h [Q,S 2 ]. 

Remark 6.3. (a) In other words, Si -< S2 if Si, S2 are local Q-droplets and Si is a local (Q, ti)-droplet 
with localization S2, where fj = f (Q, Si). 

(b) It follows from the definition that if Si, S2 are Q-droplets with Si < S2, then f(Q, Si) < f(Q, S2). 

(c) If Si < S2 and S2 < Si for two local Q-droplets, then Si c S2 and S2 c S2, and so Si = S2. 

Proposition 6.4. Let Si, S2 be two local Q-droplets with S\ c S2. Then S\ < S2 if and only ifQs l < Q 
holds on S2. 

Proof. This follows from Corollary 15. 131 After all, for local Q-droplets we do not need to check 
conditions (i)-(iv); only (v) remains. Moreover, by continuity and the fact that Q-droplets lack 
Q-shallow points, the q.e. statements hold everywhere. □ 

Lemma 6.5. For a local Q-droplet S, we have S < S. 

Proof. This follows from the definition of the "<" relation together with Proposition 15.91 □ 

There is one more property we need to check to show that "<" defines a partial ordering. 
Lemma 6.6. If Si, S2, S3 are three local Q-droplets with Si < S2 and S2 < S3, then Si < S3. 
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Proof. If we use that S 2 < S3, we see from Lemma [6.11 that SfJQ, S3] is a local Q-droplet with 

S tl (Q / S 3 )cS 2 = S f2 (Q / S 3 ). 
Using that Si < S 2 , we appeal to Lemma [572l and get 

S h [Q,S 3 ] = S h [Q r S 2 ] = S l , 

so that Si -< S 3 , as claimed. □ 

Remark 6.7. Si c S2 does not imply Si < S 2 . For example, suppose Q has two global minima at 
the points and 2, and suppose the minima are non-degenerate. Consider 

Si :=S fl (Q, Ei), S 2 :=S t2 (Q,E 2 ), where < h « t 2 « 1, 

with Ei = D(0, 1) and E 2 = D(0, 3). Then Si c S 2 but Si S 2 . This is easy to see using the 
characterization of Proposition ^. 41 

6.2. A comparison principle. We keep the setting of the previous subsection. We recall the 
definition of the polynomially convex hull phull(E) of a compact set E from Subsection l4.7l The 
set phull(E) \ E is the union of all the bounded components of C \ E. 

Proposition 6.8. Suppose Si, S 2 are two local Q-droplets with Si c S 2 . We then have Q$ 2 < Qsj on 
phull(Si), with equality on S\. Moreover, if for some z e int[phull(Si)] wehave Qs 2 (zo) = Qsi( z o)> t nen 
Qs 2 = Qsi holds on the component of int[phull(Si)] that contains Zo- 

Proof. The difference Qs 2 - Qsj is in W 2,p and therefore continuous, and it is subharmonic, as 

A[Qs 2 - QsJ = ls 2 \ Sl AQ > a.e. on C. 

Moreover, by Corollary 15.131 Qs 2 = Q on S 2 and Qsj = Q on Si, and so Qs 2 - Qsj = on Si 
as Si C S 2 . The inequality Qs 2 - Qsj < now follows from the maximum principle. The last 
assertion follows from the strong maximum principle. □ 

We see that a local Q-droplet S 2 with Si < S 2 does not grow in the direction of the interior 
holes of Si: 

Corollary 6.9. Suppose Si, S 2 are two local Q-droplets with S\ < S 2 . Wefeft/zai?e[S 2 \Si]nphull(Si) = 
and Qs 2 = Qs 1 on phull(Si). 

Proof. If Si < S 2 we have Qs t < Q on S 2 (cf . Proposition I6.4L and we also have Qs 2 = Q on 
S 2 (cf. Corollary 15.131 . In view of Proposition 16.81 it follows that if zo e S 2 n phull(Si), then 
Qs 2 (zo) = Qsi(zo). So, if z e S 2 n int[phull(Si)], another application of Proposition l6.8l shows that 
Qs 2 = Qsi holds on the component of int[phull(Si)] which contains Zo. Taking the Laplacian, we 
find that lsiAQ = ls 2 AQ a.e. on the component Comp(zo) of int[phull(Si)] which contains Zo, 
which leads to 

Si n Comp(z ) = S 2 n Comp(zo). 

Since zo e S 2 we also must have zq e Si. We conclude that [S 2 \ Si] n int[phull(Si)] = 0, and a 
fortiori [S 2 \ Si] n phull(Si) = 0. But then Q$ 2 - Qs 1 is harmonic in int[phull(Si)] and vanishes on 
<?[phull(Si)] c dS\ c Si, and the conclusion Q$ 2 - Qs 1 = on phull(Si) is immediate. □ 



COULOMB GAS ENSEMBLES AND LAPLACIAN GROWTH 



33 



6.3. Domination chains of local droplets. We are interested in chains of local Q-droplets. 

Definition 6.10. A domination chain of local Q-droplets is a (continuously indexed) family of 
Q-droplets {St}t, where the index t ranges over a nonempty interval I C R+, with left endpoint 0, 
such that t = f (Q, St) and 

h<h <^> S h < S t2 . 
The domination chain is terminating if the interval I is given by < t < t t , for some T, with 
< U < +oo, and non-terminating if it is given by < t < U for some U with < U < +oo. In case 
the domination chain is terminating, we say that it terminates at St,. 

Lemma 6.11. Given a local Q-droplet S„ there is exactly one domination chain of local Q-droplets that 
terminates at S*. 

Proof. By Lemma l6Tl St := St[Q, S,] for < t < f, := t(Q, S») defines a continuously indexed 
collection of local Q-droplets, and by Lemma [5.2| it is a (terminating) domination chain. Finally, 
if Sjj is a local Q-droplet with S» < S,, then by definition, it is of the form Sj = St„[Q, S*] with 
:= KQt ^ U, so the domination chain is unique. □ 

6.4. Maximal domination chains of local Q-droplets. We keep the setting of the previous sub- 
section. We shall need the concept of a maximal domination chain of local Q-droplets. 

Definition 6.12. A domination chain of Q-droplets is maximal if it is contained in no larger 
domination chain of local Q-droplets. 

Maximal domination chains of Q-droplets can be either terminating or non-terminating. If the 
chain is indexed by the unbounded interval I = R+ then it is automatically non-terminating. If 
the chain is indexed by a bounded interval, then it can be non-terminating only if the droplets 
develop "arms" or "islands" that tend to infinity: 

Theorem 6.13. Let [S t }m be a maximal non-terminating domination chain of local Q-droplets. Then the 
union 

Su :=|J St 

is an unbounded subset ofC. 

Proof. We suppose Su is bounded, and form S, = clos Su, which is then compact. We are to show 
that the non-terminating domination chain {St}tei cannot be maximal. The interval J is given by 
< t < U for some t» with < t, < +oo. For tel, we let at be the positive measure dot - Is, AQdA, 
which has total mass \\o t \\ = t. Let a, be given by da, = ls u AQdA, which has total mass ||a,|| = t». 
Then Ot — > a, in norm as t — > f„ and in fact the corresponding densities converge in LP: 

l s , AQdA -> l Su AQdA in U(C) as t -> U. 

By the well-known properties of the 2D Hilbert transform, we find that the associated potentials 
converge in W 2 ' p : U^' Sl — > U®' Su as t — > t,. Also, we easily check that if the constants y*{Q, S t ) and 
y*(Q, Su) are as in 1 15.31 1, we have y*(Q, St) — > ~y*(Q, Su) as t — > £,. As a consequence, 

Q St = f t {Q, St) - !iQ' s ' -» Q Su = yl(Q, S t ) - U& s " in W 2 * as t -» £,. 

By Sobolev imbedding the convergence is locally uniform. Since Qs, = Q on St we get that 
Qsu = Q on Su- By continuity, then, we find that Qs u = Q on S, = clos Su- Next, by [18J, p. 53, 
we see that AQs v = AQ a.e. on S,, that is, ls u AQ = AQ a.e. on S„. Expressed differently, we have 
lg u AQ = ls.AQ as elements of LF(C). In particular, AQ > holds a.e. on S». By construction, S, 
has no Q-shallow points, a property this set inherits from the individual droplets St, t e I. In view 
of Corollary l5.17[ S, is a local Q-droplet. It remains to show that we may add S, as a terminal local 
Q-droplet for the domination chain, thereby defeating the maximality of the non-terminating 
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domination chain. To this end, it suffices to obtain that S t < S, for tel. We pick a f with 
t < f < U, and use S t < S? to deduce that Q$ t < Q on S? (Proposition 16 ,4|l . By letting V — > t», we 
get that Qs t < Q on Su, and by continuity that Qs t < Q on S„. By Proposition 16.41 this means that 
Sf < S*. The proof is finished. □ 

The following definition is useful. 

Definition 6.14. A local Q-droplet S is maximal if for any other local Q-droplet S' the relation 
S<S' implies that S = S'. 

Corollary 6.15. A maximal domination chain {S t }tei of local Q-droplets either terminates at a maximal 
local Q-droplet, or is non-terminating, in which case the set S u ofTheorem \6.13\ is unbounded. 

6.5. Richardson's formula. We keep the setting of the previous two subsections. We would like 
to understand the flow evolution t i-» Sf of domination chains (or containment chains, see the 
next subsection) of local Q-droplets. A natural way to do this is to analyze the effect of the flow 
when we use harmonic functions as test function (i.e., we calculate "harmonic moments"). 

Proposition 6.16. Suppose S,S' are two local Q-droplets with S c S'. Then for all h e W 2/1 (C) (local 
Sobolev class) that are harmonic in C \ S and bounded near infinity, we have (with t = t(Q,S) and 
t' = KQ,S')) 

1 hAQdA = (t' -f)/z(oo). 



f ■ 

JS'\S 



Proof. The formula holds for constant h, by the choice of t, t'. So, by subtracting a constant, we 
may take h(oo) = 0. We have 

(6.1) f JiAQdA = fh l s AQdA = f hAQ s dA = f Q s AhdA, 

Js Jc Jc Jc 

and the analogous identity holds for S' as well. By forming the difference between | |6.1t for S and 

S' we see that 



f hAQdA = f \Q s ,-Q s ]AhdA = 0, 
JS'\S Jc 



because Ah = on C \ S while Qs' -Qs = Q- Q = 0onS (see, e.g., Corollary |5.13t . To finish the 
proof, we just need to justify f6.lt . By Green's formula, we have 

(6.2) f \hAQs - Q s Ah]dA = 2 f \hd n Q s - Q s d n h]ds, 

Jd(o,r) Jt(o,r) 

where ds is normalized arc length (i.e., arc length divided by 2n) and d n is the exterior normal 

derivative. Next, we observe that as |z| — » +oo, we have the asymptotics 

h = 0{\z\~ 1 ), |V/i| = 0(|z|- 2 ), Q s =0(log|z|), |VQ S | = Odzp 1 ), 

because both h and Qs are harmonic in C \ S with given asymptotical behavior. By letting _R — » +oo 
in | |6.2t , we obtain f6.lt . The proof is complete. □ 

Corollary 6.17. Suppose S, S' are two local Q-droplets with S c S'. Also suppose that the interior int S 
has finitely many components. Then for all h continuous and bounded in C \ int S, which are harmonic in 
C \ clos int S, we have (with t = t(Q, S) and f = t(Q, S')) 

\ hAQdA = (f -t)h(oo). 

JS'\S 

Proof. By Mergelyan-type approximation we can find a sequence of bounded C°°-smooth func- 
tions h n that are are harmonic in C \ clos int S, such that h„ — > h uniformly on C \ int S as n — > +oo. 
The assertion now follows from Proposition 16.161 □ 
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6.6. Differential form of Richardson's formula. We keep the setting of the previous subsections, 
and introduce the concept of a containment chain. 

Definition 6.18. A containment chain of local Q-droplets is a (continuously indexed) family of 
Q-droplets {St)t, where the index t ranges over a nonempty interval / c R+, with left endpoint 0, 
such that t = t(Q, S t ) and 

h<t 2 <=> S h c S f2 . 

Let {Sf}f e ; be a containment chain of local Q-droplets, and let U denote the right endpoint of I 
Let I~ be the interval obtained from / by removal of U (if U & I, we put I~ :- I). 

Lemma 6.19. The map 1 1-» St, t e I, is continuous in the Hausdorff metric except for a countable subset 
of the interval I. 

Proof. For to e I~, we form the compact sets 

s-=cioslJs t , s; o = f)s t . 

t:t<t t:t>t 

Then S~ c S t c S+ holds for each for t e 1~. Note that S~ is well-defined also when to = U. For 
t e I~, we put 

5~(t) := maxdistc(z, SJ), 5 + (t) := maxdistc(z, S f ). 

zeS, zeS+ 

Next, for positive e, we consider the sets 

DJ := {ter: 6~(t) > e}, := (fef: 6 + (f) > e}. 

We argue that for each positive e, the sets D~,D* are countable and that the only possible 
accumulation point is U (and if U is an accumulation point, then the set Sj must be unbounded). 
Indeed, if, for some f with < f < U, the set D~n]0, t'] has N = N(t',e) elements, then the local 
Q-droplet Sf contains N points which are e-separated (the distance between any two different 
points is at least e). We get an effective bound on N in terms of the diameter of the compact set 
Sf. The analogous argument applies to in place of D~. a 

We let gj^ denote harmonic measure for the open set C \ S t with respect to the point at infinity. 
This is a probability measure whose support is contained in <9phull(Sf) C dSt (the effect on a test 
function is that we get the value at infinity of the harmonic extension). 

Proposition 6.20. Suppose the map 1 1-> St, t £ I, is right continuous at to 6 I~. Suppose moreover that 
int Sf has finitely many components and that St = clos int Sf . Then for all g e C(C) we have 



lim -L- f g AQdA = f gdcv 
t->q t - to J SA s J ds 



(to) 

oo / 



that is, we have the weak-star convergence of measures 



lim -±-l SASlo AQdA = dJi°\ 
t->t+ t — to " 



Proof. Let h denote the function which coincides with g on Sf and extends harmonically (and 
boundedly) to C \ Sf, so that in particular 




Then h is continuous and bounded in C (see, e.g. [ 10 1 for a discussion of the Dirichlet problem). 
Since S fo = clos int S fo , Corollary 16 . 1 71 applied to to h gives 

-L- f gAQdA = h(oo) + -±- f (g-h)AQdA. 
t - f o Js t \s, t - to Js,\s t0 
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It remains to show that the last term on the right hand side tends to zero as t — > to- This 
follows from the fact that h(z) - g(z) — > as z — > dSt and that St \ St by the right continuity 
assumption. □ 

Remark 6.21. Proposition ^. 20| states that (under regularity assumptions) the infinitesimal growth 
of the local Q-droplets is in the exterior direction only. If the containment chain of local Q-droplets 
were to grow in the direction of the internal holes, the containment chain could not possibly be a 
domination chain (cf. Corollary 16.91 . 

6.7. Richardson's inequality. We now show that under modest regularity conditions, contain- 
ment chains of local Q-droplets are in fact domination chains. 

Theorem 6.22. Suppose S, S' are two local Q-droplets, with S c S'. Then the following are equivalent: 

(i) S < S'. 

(ii) For all real-valued functions h e W 2,1 (C) (local Sobolev class) that are subharmonic in C\S, harmonic 
in C \ S', and bounded near infinity, we have (with t = t(Q, S) and t' = t(Q, S')) 

(t'-t)h(oo)< ( hAQdA. 

JS'\S 



Proof. We first show that (i) => (ii). We note that the inequality is an equality when h is constant 
(see, e.g., Proposition 16 . 1 6) . This allows us to restrict our attention to h with /?(oo) = 0. As in the 
proof of Richardson's formula (Proposition 16 . 1 6t , we find that 

f hAQdA = f [Qs- - Qs]AhdA = f [Q s , - Q s ]AhdA > 0, 

JS'\S JC JS'\S 

since Qs < Q = Qs', by Proposition 16 .41 and Corollary 15. 131 

We turn to the implication (ii) => (i). We take h(oo) = 0, and get (as above) from (ii) that 

(6.3) 0< f hAQdA = f [Q s , - Q s ]AhdA. 

JS'\S JS'\S 

The question now is what kind of functions Ah are possible here. We have automatically Ah e 
L X (S') while Ah > a.e. on S' \ S. We also need to impose that 



X 



AhdA = 0, 

IS' 

as a consequence of the behavior of h near infinity. In fact, any real-valued function g e D(S') for 
some q with 1 < q < +oo with g > a.e. on S' \ S and 

(6.4) J gdA = 0, 

is of the form g - Ah for an h as in (ii) with h(oo) = 0. It follows from | |6.3t that 

(6.5) 0< f [Q S '-Qs]gdA. 

JS'\S 

As it is easy to fulfill | |6.4t by placing an L'-integrable negative mass on S to compensate for the 
positive mass on S' \ S, on S' \ S the function g is basically any positive 13 function on S' \ S. This 
is only possible if Q s < Q s > on S' \ S, and as Q s - = Q on S', we get Q s < Q on S' \ S. Since Qs < Q 
holds automatically on S, we see that Qs < Q on S'. The conclusion S < S' now follows from 
Proposition [6^4] The proof is complete. □ 

The proof of Theorem |6.22| has the following consequence. 

Corollary 6.23. Suppose S' is a local Q-droplet, and that S c S', where S is compact and lacks Q-shallow 
points. If condition (ii) of Theorem \6.22\ is fulfilled, then S is a local Q-droplet, and S < S' . 
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Proof. As in the proof of Theorem |6.221 we get from condition (ii) of that theorem 
< f hAQdA = f [U Q ' S - U Q ' s ']AhdA = f [U as - U Q ' s ']AhdA + \ [U Q ' S - U Q ' s ']AhdA, 

JS'\S JC JS JS'\S 

provided that h e W 2/1 (C) (local Sobolev class) is subharmonic in C \ S and harmonic in C \ S', and 
bounded near infinity, with h(co) = 0. As in the proof of Theorem 16.221 we choose h as (minus) 
the logarithmic potential of g, where g e L''(S') has g > on S' \ S and 



X 



gdA = 0. 

Js> 
so that 

(6.6) < f [U Q ' S - U Q ' s ']gdA + f [U as - U Q ' s ']gdA, 

JS JS'\S 

If we choose g such that g = on S' \ S, we have equality (since then the inequality applies to —g 
as well): 



[ U Q,s _ u as ']gdA = 0. 

>s 

As g is now arbitrary except that its integral over S vanishes, we conclude that U® ,s - U^' s ' is 
constant on S. Call the constant c: U®' S = c + JJ Q S on S. Since S' is a local Q-droplet, the function 
jjQ'S' + Q is constant on S' (cf. Corollary 15.171 and consequently, U^ ,s is constant on S (after all, 
S c S'). We conclude that S is a local Q-droplet. That S < S' now follows from Theorem |6.22| □ 

Theorem 6.24. Let {S t }tei be a containment chain of local Q-droplets. If, for almost every t e I, the set 
int St has finitely many components and St = clos int St, then {St}t€i is a domination chain. 

Proof. We consider the W 2,p -smooth function 



V{Z,n;t):=U^ s >{E,)-U^ s >( n )= f log 

Js, 



z - rj 



AQ(z)dA(z), t 6 I, 



and note that 

(6.7) V(^rj;t) = Q St (il)-QsM), tel. 

Let \i, v be two compactly supported Borel probability measures which are absolutely continuous 
with densities in U for some q, 1 < q < +oo. We need the expression 



V^(t):= f V{^,q;t)d^)dv{q)= f U v ^AQdA ; f 6 J, 
Jc 2 Js, 

where !i v_f ' := IF - U 1 ', and U'', !J V are the usual logarithmic potentials. The functions W l , IF are 
in W 2 <? and therefore continuous (and bounded). The function !i v_f ' is harmonic off supp(v - /i), 
and its value at infinity is IF - ^ 00 ) = 0. We have 

(6.8) V' v (f) - V^' v (t) = f LT^'AQdA, for t,t' e I with f < f, 

Js,,\s, 

which gives 

W^ v (f) - V^'(t)\ < IILf^ll^^) f AQdA = (£' - t)ll U l '-q L ~ (C) , for t, f e J with f < £', 

Js,,\s, 

since AQ > a.e. on a local Q-droplet. It follows that the function W' v is Lipschitz continuous, 
and therefore differentiable almost everywhere. In view of ( 16.8b , its right derivative is 

[V» v ]'(t + ) = lim -L- f LT-f AQdA = f !J v -f'da;2, 
t'^t + f - f J SASi J dST 
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by Proposition 16.201 with the possible exception of a countable set of f's. If now supp v c Sf, the 
function U v ~^ 1 becomes subharmonic (and bounded) in C \ Sf, so by the maximum principle 

= LT-''(co) < f W^dco®. 

JdS™ 

We conclude that [V f, ' v ]'(f) > for a.e. t with supp v c S t . Put 

t v :- infjf e I : suppv c St], 

and note that for t e / with t > t v we have [V i/V ]'(i) > almost everywhere, and hence V f,,v is 
increasing on that sub-interval: 

(6.9) V^' v (t) < V^' v {f) for f,f 6 I with t v < t < f. 

Next, we let the probability measures \i, v get more and more concentrated, so that supp — » {£} 
and supp v — » {r|}. The inequality l|6.9t survives the limit process, and we obtain that 

(6.10) V(£, rj; t) < V(£, rj; t') for t, t' e I with t ( < t < t' , 
where 

t s := infjf e I : £ e S f }. 

The short argument which justifies this involves choosing the support of v cleverly, and this is 
made possible by the fact that a local Q-droplet lacks Q-shallow points. If we use \6.7\ , we see 
that (16. 10b expresses that 

(6.11) QsM) - QsAZ) ^ QsAn) - QsAO for t,t'ei with t 4 < t < r. 

Since for ti < t < f we have E, e St c Sf , we get that (cf. Proposition ^. 5t 

Q S( (£) = Q S( ,(£) = Q(£), 

so that l|6.11t simplifies: 

Qs t (r?) < Qs,, (rj) for t,t' e I with f { < f < f . 
By making clever choices of the point £ we can get t{ to be as close to as we need, and so 

Q SI (V) < QsA 1 !) for t,t' el with f < r. 
For T] e Sf we have Qs t , (/]) = Q(f]), and we derive that for t,t' e I with f < t', we have 

Qs,('7) < Qfa) r?eSf. 
By Proposition ^. 41 we get Sf < Sf for all t,t' e I with f < f , and {S{}f S j is a domination chain. □ 

7. The Hele-Shaw equation 

7.1. Smooth curve families (laminations). We need the following definition. 

Definition 7.1. A family of simple curves Tt (where t runs over some interval) in C is a C°°-smooth 
lamination if 

(i) T f n I> = holds for t + f, and 

(ii) Each curve has a local parametrization z = yt(6) (6 runs over some interval), such that the 
function y(6, t) :- yt(6) is a local C°°-diffeomorphism. 

We will alternatively use the term C^-smooth curve family as synonymous to C°°-smooth lami- 
nation. We mention that it is of course also possible to define laminations with a lower degree of 
smoothness than C°°. The normal velocity v„ = V„(z), z 6 T t , may be defined as follows: 

1 

v„ := (d t y,n) = — — ■ lm[d t yd g y], 
\d e y\ 
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where the inner product is that of C = R 2 and n is a unit normal to T t . It is easy to see that the 
definition does not depend on the choice of parametrization y. Indeed, if we write 

f t (9) = y(9,t) := y(9(9,t),t), 
where 9 i-» 8(9, t) is a local diffeomorphism, then 

d t y = dey d t 6 + d t y, d$y = dey d$d, 

so that 

1 1 

-3— lm[d t yd 9 y] = — — 5— Im \\d g y\ 2 d t 6d s d + d t yd g fd s e] 
\d$y\ \d g yd s e\ 

= T3^HT 73— r IMdtydef] = ±7T— r hn[d t yd e f], 
\dsd\ \dey\ \d e y\ 

where there is a sign change if the coordinate change reverses the direction of the unit normal 
vector. 

Lemma 7.2. Let T t be a C^-smooth lamination of Jordan curves, such that the domain D t interior to T t 
increases with t. Then, for continuous f : C — » C, we have 

d_ 

dt 



\ fdA = 2 f fv n ds. 



Proof. We identify the area form with the area measure according to, e.g., dz A dz = 27iidA(z). We 
may assume that for t, to close to one another with to < t,D t \ D t(l is parametrized by y T (6) = y(0, t) 
where < 6 < 1 and to < % < t, with periodicity boundary conditions in 6: y(0, t) = y(l, t). We 
let R(t , t) denote the rectangle [0, 1] x [to, t], so that 



f fdA=^- [ f(y(d,T))d 7 Ady. 

JD t \D t „ Zm JR(t ,t) 



We calculate: 

dy A dy = [deyd t y - dgyd t y] dd A dt = 2ilm[deyd t y] dd A dt = 2i\d0y\v n dddt, 

where we have identified a form with the corresponding measure. We identify |(?ey|d(9 as arc 
length along T t , so that |(?ey|d(9 = 27ids(0), and therefore, 

1 



_ :dy A dy = 2v n ds(d)dt. 
2tzi 



The assertion is now immediate. 



7.2. The Hele-Shaw flow equation. We assume we have a C°°-smooth lamination of Jordan 
curves Tt, and let Dt denote the interior domain while Qf is the exterior (unbounded) domain. 
We also write K t := clos D f = C \ Q f , so that K t is compact. The classical Hele-Shaw equation relates 
the normal velocity v n to the normal derivative of the Green function (for the Laplacian) of the 
exterior domain Q f when one of the two coordinates is the point at infinity (the factor | comes 
from our choice of normalizations): 

(7.1) v n = \d n G t on T t , where G f = G(-,oo;Q f ). 

The Green function Gf is always positive in Q f and vanishes along the boundary Tt, and n is 
taken in the exterior direction, so that d n Gt > on Tt- Actually, d n Gt is the Poisson kernel of Q 
for the point at infinity, so that hd n Gt times normalized arc length measure has the interpretation 
of d<jj^, harmonic measure at infinity for the domain O t . There is also a weighted analog of 1)7. lb : 
the weighted Hele-Shaw equation is 

(7.2) pv n = \d n G t on r f . 
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The function p is the weight, and it is assumed to be C°°-smooth with p > point-wise. It is 
possible to interpret the introduction of the weight as a change of the geometry (cf. Ifl6l , [14J, 
Ifl5l ). In the sequel, we will use p = AQ. 



Definition 7.3. We say that an increasing family of compact sets {K t } t (where t ranges over some 
interval) is a generalized solution of the weighted Hele-Shaw equation with weight p - AQ if 

(i) AQ > on U t K t , if 

(ii) for each t e I, K t lacks Q-shallow points, and if, 

(iii) for all / 6 C(C), the function 



t i-» 



f /AQcM 

Jk, 



is absolutely continuous and for a.e. t we have (Tt = dQf where := C \ phull(X f ) is the 
unbounded component of the complement C \ Kt and cu^, is harmonic measure at infinity for ) 



^ J fAQdA = £ fdv®. 



Note that no smoothness requirement is imposed on the compact sets Kt as in the standard 
formulation of the weighted Hele-Shaw equation l|7.2t . The way things are set up, strong solu- 
tions of the weighted Hele-Shaw equation (i.e., solutions of j7.2\ ) are automatically generalized 
solutions. In short, the equation asks that the compact sets Kt grow according to the law 

^[l Kl AQdA] = dJl 

Proposition 7.4. Let {K t }t be an increasing family of compact sets, where t ranges over an open interval 
I, and suppose AQ > on D te jKt, and that K t lacks Q-shallow points, for each tel. Then {K t } t<E i is a 
generalized solution of the weighted Hele-Shaw equation with weight AQ if and only if, for all t,V el 
with t < t' , and for all real-valued f e C(C), 



f fAQdA = f f fdco^dx. 

JK t ,\K, Jt JT, 



Proof. This is just an application of Calculus. □ 
So, the weighted Hele-Shaw equation corresponds to the disintegration of measures 

l KAKl AQdA = £ do^dx. 

It follows from the standard properties of the harmonic measure that if / e C(C) is real-valued, 
bounded, and subharmonic in C \ K t while it is harmonic near infinity, then 

(7.3) f fAQdA = f f fdJZ ] dT> f f(oo)dz = (f - f)/(oo). 

JKf,\K, Jt Je t Jt 

This strongly resembles Richardson's inequality for local Q-droplets (Theorem 16.221 . The com- 
parison with Theorem |6.22| suggests the concept of a weak solution to the Hele-Shaw equation. 

Definition 7.5. We say that an increasing family of compact sets {Kt}t (where t ranges over some 
interval) is a weak solution of the weighted Hele-Shaw equation with weight AQ if 

(i) AQ > on U t K t , if 

(ii) for each t e I, K t lacks Q-shallow points, and if, 

(iii) for all real-valued / e W 2,1 (C) (local Sobolev class), 

(f - f)/(oo) < I fAQdA for t,t' el with t < t' , 

Jkak, 

provided / is subharmonic in C \ K t , harmonic in C \ Kt>, and bounded near infinity. 
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Proposition 7.6. A generalized solution of the Hele-Shaw equation is a weak solution. 

Proof. It is known that it suffices to have condition (iii) of Definition 17.51 fulfilled for / e W 2,< ? for 
some q slightly bigger than 1 . Such functions are continuous, so the assertion is immediate from 
ES. □ 



We note that the sets K t need not be local Q-droplets, although that is one particular instance. 
The analogy with that case suggest the following. 

Definition 7.7. An increasing family of compact sets {K t } t is correctly indexed if 

t = t(Q,Kt)= f AQdA. 
Jk, 



This is in agreement with (|7.3t (or with Definition 17.51 for / = 1 (since the inequality applies 
to / = -1 as well the inequality is of course an equality). 

It is known that the Hele-Shaw equation behaves like the heat equation, in that one direction of 
time t is stable and the other is unstable. Here, the evolution f i-» Kt is unstable when t increases, 
and stable when t decreases. 

Theorem 7.8. Let K, c C be compact, with AQ > a.e. on K t . We assume that U ■= t(Q,K t ) > 0, and 
that K, lacks Q-shallow points. Then there exists a correctly indexed weak solution 1 1-» K t of the weighted 
Hele-Shaw equation with weight AQ on the interval < t < U, such that K t , = K t . The solution is unique. 

Proof. We consider the function Q, := -U®' K, r where U® ,K ' is as in | |5.2t . It has 
AQ, = Ik, AQ, Q,(z) = t, log |z| 2 + O(l) as \z\ -» +oo, 

and we can define 

(7.4) K t :=S t [Q„K,], 0<t<t„ 

and K tt := K,. The way things are set up, K> becomes a Q» -droplet (cf. Corollary |5.17t . Moreover, 
in view of Lemma 16.11 (with Q» in place of Q), the sets Kt are local Q» -droplets. We see from 
Theorem 16.221 that the sets K t form a domination chain of local Q» -droplets if and only they form 
a weak solution of the Hele-Shaw equation. 

It remains to establish that the weak solution t t-^> K t unique. So, suppose t i-» Kt is a weak 
solution, which need not be of the form (17.41 . We claim that Kt is a Q» -droplet for < t < f». We 
know that K, is a local Q» -droplet, that Kt has no Q-shallow points, and that Kt c K,. In addition, 
the weak solution condition entails 

(f„ - f)/(oo) < f /AQdA, 0<t<U, 

Jk,\k, 

provided / e W 21 (C) (local Sobolev class) is real-valued, subharmonic in C \ Kt, harmonic in 
C \ K,, and bounded near infinity. An application of Corollary 16.231 shows that K t must also be 
a local Q» -droplet, with Kt < K, with respect to the weight Q». The uniqueness part is now a 
consequence of Lemma |6.11| □ 

Remark 7.9. (a) A key element of the proof of Theorem 17.81 is the identification of the weak 
solutions of the Hele-Shaw equation t i-> Kt with domination chains with respect to the weight 

Q- 

(b) Theorem 17.81 supplies existence and uniqueness in the backward time direction. It is not 
difficult to see that there is even local uniqueness in the backward time direction. However, in 
the forward time direction, there is generally neither existence nor uniqueness. An example of 
non-uniqueness can be based on, e.g., the setting of Remark |6. 71 We now discuss non-existence. 
In the context of Theorem |7.81 the difference U®' K " - U® ,Kt is constant on Kt for < t < t, (see, e.g., 
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the proof of Corollary 16. 23t ; let c(t) be that constant. We consider the functions H t :- Q - U^ ,Kt 
and H, := Q - U®' K ", which have AH t = a.e. on K t and AH, = a.e. on JC, respectively. We 
have H f - H» = U Q - K - - U Q ' K > = c{t) on K t , we we write as H t = c{t) + H, on K t . So H< restricted 
to K t is supposed to have an extension to K, - the function H t := c(t) + H, - with AHj = a.e. on 
K*. This adds an additional smoothness requirement on H t for < t < U, which suggests that 
Kt cannot be an arbitrary compact subset of C with AQ > a.e. on Kt which lacks Q-shallow 
points. But K t is uniquely given for < t < U (the backward direction) for arbitrary compacts K, 
lacking Q-shallow points. So with very irregular K, we should be able to arrange that we have 
non-existence in the forward time direction. Another reason for non-existence in the forward 
direction is the existence of maximal local Q-droplets (see the next section for details), at least for 
some Q with AQ = 1. 

A proof of the following statement can be based on Proposition ^. 201 The only part that needs 
checking is the absolute continuity requirement, which we leave to the interested reader. 

Proposition 7.10. Suppose t 1-> K t is a weak solution to the Hele-Shaw equation and that for almost all t 
the sets int K t have finitely many components. Then t^>K t is a generalized solution. 
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